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ABSTRACT 

This  report  covers  analytical  and  experimental  studies  of  the  behavior  of 
burning  Composition  B  and  PBX  9404  under  various  degrees  of  confinement. 
These  studies  were  undertaken  to  determine  the  properties  of  the  gas,  the 
conductive  heat  transfer  to  the  explosive,  the  radiant  heat  transfer  to  the 
surfaces  of  the  high-explosive  and  chamber,  and  the  heat  evolved  by  the 
burning  of  the  explosive.  The  work  resulted  in  a  theory  which  describes  the 
burning  of  the  explosives  under  unconfined  and  confined  conditions,  and  the 
ensuing  pressure  buildup  during  the  period  the  explosive  remains  intact  and 
the  reaction  is  limited  to  its  surface.  The  results  are  in  good  agreement 
with  experimental  measurements  of  pressure,  gas  temperature,  chamber 
temperature,  and  heat  flux  to  meters  in  the  gas  stream.  Beyond  a  pressure 
of  about  1000  psia,  the  samples  of  Composition  B  cracked  causing  an  abrupt 
increase  in  the  rate  of  pressure  rise.  The  cause  of  these  phenomena  and 
the  question  of  whether  or  not  PBX  9404  cracks  were  not  resolved.  The 
results  of  this  work  may  be  used  to  predict  the  minimum  order  of  explosion 
which  may  arise  from  a  system  that  is  exposed  to  heating. 
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1.  INTRODUCTION 

Numerous  accidents  have  occurred  in  which  nuclear  weapon 
systems  have  become  involved  in  fire.  From  1956  to  1963  a  number 
of  weapon  systems  have  been  purposely  exposed  to  fire  to  ascertain 
the  behavior  of  the  high  explosive  and  the  electronic  circuits.  The 
results  of  these  tests  are  presented  in  the  DASA  Reports  1060.  1066, 
1234,  and  1276  (see  Bibliography),  For  purposes  of  safety,  infor¬ 
mation  is  generally  required  pertaining  to  the  minimum  ignition  time 
of  the  high  explosives,  the  likelihood  of  a  high-explosive  detonation, 
and  the  likelihood  of  a  nuclear  detonation. 

The  ignition  times  of  nuclear  weapon  systems  have  been  pre¬ 
dicted  to  within  20  percent.  This  was  accomplished  Ly  building 
models  of  critical  regions  of  the  weapon  system  and  exposing  the 
models  to  the  desired  rate  of  heating.  A  detailed  discussion  of 
modeling  techniques  is  presented  in  the  DASA  Report  1130. 

The  likelihood  of  a  nuclear  detonation  from  heating  of  the 
high-explosive  is  dependent  on  the  behavior  of  the  various  electronic 
circuits  and  on  the  condition  of  the  HE  system  at  the  time  of  delivery 
of  any  critical  electrical  signal.  DASA  Report  1276  describes  the 
behavior  of  a  number  of  electronic  circuits  within  actual  weapon 
systems  that  were  exposed  to  fi-e.  Detailed  laboratory  studies  of 
the  hehavior  of  the  various  electronic  components  in  high  temperature 
environments  have  been  conducted  by  organizations  such  as  Sandia 
Corp.  .  Albuquerque.  New  Mexico.  From  such  work,  it  is  possible 
to  conservatively  estimate  the  probability  of  a  nuclear  detonation  of 
a  given  weapon  system  in  a  fire  environment. 

Work  has  also  been  conducted  to  develop  means  of  predicting 
whether  a  given  weapon  system  would  or  would  not  explode  in  a  given 
fire  environment.  This  work  was  initiated  by  a  feasibility  study  in 
1956  by  I  I  T  Research  Institute  to  explore  the  possibility  of 
developing  means  of  predicting  HE  explosions  and  means  of  pre¬ 
venting  such  explosions.  This  study  indicated  that  such  an  endeavor 
appeared  promising  and  led  to  the  studies  covered  in  this  report. 

The  results  of  the  initial  study  are  presented  in  the  DASA 
Report  1  234.  In  this  study,  cylinders  of  Composition  B  were  burned 
at  the  bottom  of  tubular  chambers.  Vents  were  provided  at  the  top 
of  the  chamber  to  vary  the  rate  of  pressure  build-up.  This  work 
revealed  that  the  behavior  of  burning  explosives  under  pressures  of 
1  to  2  atmospheres  is  reasonably  reproducible  and  therefore  amen¬ 
able  to  prediction.  Measurements  of  thermal  irradiation  at  various 
locations  within  the  gas  stream  indicated  that  thermal  radiation  is 
the  principal  source  of  heat  flow  to  the  thin  layer,  i.  e.  ,  approxi¬ 
mately  0.  05  inches,  of  explosive  undergoing  reaction.  Thermal 
radiation  is  of  less  significance  at  higher  pressures  of  the  order 
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cf  tens  of  atmospheres  due  to  the  rapid  increase  in  the  heat  conduction 
from  the  hot  gas  with  pressure.  It  was  found  that  vent  holes  having 
a  cross-sectional  area  of  about  1  percent  or  more  of  the  area  of  the 
burning  explosive  were  sufficient  to  maintain  near-atmospheric  pres¬ 
sures  in  the  chamber. 

This  report  covers  the  experimental  and  theoretical  studies 
of  the  behavior  of  burning  Composition  B  and  PBX  9404  under  pres¬ 
sures  as  high  as  several  hundred  atmospheres.  Section  2  presents 
the  significant  properties  of  the  explosive  and  test  chamber  required 
for  the  theoretical  studies  and  develops  a  theory  describing  the  flow 
of  thermal  radiation  from  the  hot  gas  stream.  Section  3  discusses 
the  experimental  burning  of  the  high  explosives  under  pressures  of 
1  to  3  atmospheres.  In  addition,  several  analyses  were  conducted 
to  predict  the  burning  velocity,  the  gas  temperatures,  the  thermal 
irradiation  from  the  gas  stream,  the  temperatures  of  the  test  chamber 
and  the  gas  flow  through  the  vent  holes.  A  comparison  of  the  results 
of  the  experimental  and  analytical  studies  ranges  from  good  to  excel¬ 
lent. 


Section  4  discusses  the  burning  of  Composition  R  and  PBX  9404 
under  pressures  as  high  as  several  hundred  atmospheres.  This  was 
achieved  by  burning  the  explosive  material  in  closed  chambers  which 
were  held  together  by  rods  of  varying  strength.  The  void  volume 
within  the  chamber  was  varied  to  explore  its  effect  on  the  pressure 
history  of  the  gas.  The  results  of  the  analytical  studies  covered  in 
Section  3  were  used  to  predict  the  pressure  history  of  the  gas  evolved 
by  burning  Composition  B  and  was  found  to  be  in  good  agreement 
with  the  experimental  data  for  pressures  as  high  as  about  1000  psia. 

At  pressures  in  excess  of  this  level,  the  rate  of  pressure  rise  in¬ 
creases  very  abruptly.  The  cause  of  this  abrupt  pressure  rise  is 
attributed  to  the  increase  of  burning  area  caused  by  cracking  of  the 
Composition  B.  Insufficient  data  were  found  to  conduct  a  similar 
analysis  of  the  burning  of  PBX  94  04. 
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2.  BASIC  DATA  AND  FORMULAS  PERTAINING  TO  THE  BEHAVIOR 

OF  BURNING  COMPOSITION  B  AND  PBX  9404 

This  section  covers  the  basic  data  and  formulas  required  for 
subsequent  theoretical  studies  which  describe  the  behavior  and  effects 
of  the  burning  explosives  under  various  environments.  The  theoretical 
studies  will  be  presented  later  in  the  report. 

a.  Heat  Released  by  Burning  Composition  B  and  PBX  9404 

The  heat  released  by  the  burning  of  a  given  mass  of  explosive 
plays  a  sensitive  role  in  affecting  the  heating  and  rate  of  burning  of 
high-explosive  materials.  A  preliminary  study  of  the  heat  released 
by  burning  Composition  B  under  pressures  of  1  to  2  atmospheres  is 
presented  in  the  DASA  Report  1234.  Additional  data  were  acquired 
under  the  present  program  from  a  literature  study  and  by  experi¬ 
mentation.  In  these  experiments,  a  given  mass  of  Composition  B 
or  PBX  9404  was  burned  to  completion  in  the  tubular  chamber  shown 
in  Figure  1.  The  chamber  was  insulated  with  a  light-weight  insulation 
( rockwool)  ,  to  deter  the  escape  of  heat  from  the  chamber.  The  top 
of  the  chamber  was  partially  closed  by  a  disk  that  was  formed  from 
high-temperature  insulating  brick  to  minimize  radiant  heat  loss. 

A  large  hole  was  provided  through  the  center  of  the  disk  to  allow  the 
gases  to  escape  at  a  sufficiently  fast  rate  so  that  a  near-atmospheric 
pressure  was  maintained  within  the  chamber.  The  explosives  were 
ignited  by  preheating  the  exposed  surface  with  a  122 -watt  electrical 
heater  for  15  minutes  and  passing  a  current  of  10  amperes  through 
a  0.  56  ohm  wire  embedded  across  the  face  of  the  explosive.  The 
burning  velocities  of  the  Composition  B  and  the  PBX  9404  were  2.  5 
and  5.  0  ft/hr  respectively,  and  were  found  to  be  consistent  and  uni¬ 
form  throughout  the  tests.  Most  of  the  heat  released  by  the  burning 
was  either  absorbed  by  the  steel  chamber  or  escaped  with  the  gas. 

Small  corrections  in  the  data  were  made  for  the  heat  absorbed  by  the 
insulating  materials,  for  the  heat  generated  by  the  electric  heater 
and  igniter,  and  for  the  energy  associated  with  the  velocity  of  the 
escaping  gas. 

The  heaf  absorbed  by  the  chamber  from  the  burning  of  each 
pound  of  high-explosive  was  determined  by  measuring  the  temper¬ 
atures  and  specific  heat  of  the  chamber  and  substituting  these  data 
into  the  following  formula: 

Q(  chamber)  =  i  \  C,„  ( T  -  TJ  dv  (1) 

o 
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where 

Q  (chamber) 

M 

C 


v 


heat  absorbed  by  steel  chamber  from  the 
burning  of  a  pound  of  explosive>  B/lb. 
total  mass  of  explosive,  lb 
specific  heat  of  steel  =  0.  12  B/lb*R 
temperature  of  an  increment  of  the  chamber, 
*R 

initial  temperature  of  the  chamber,  *R 
density  of  steel  =  491  lb/ ft* 
volume  of  an  increment  of  the  chamber,  ft3 
total  volume  of  chamber,  ft3 


A  similar  calculation  was  made  to  determine  the  heat  absorbed 
by  (lie  insulating  materials  Q  (insulation)  using  their  temperature 
rises,  specific  heat,  densities,  and  volumes.  The  heat  which  escaped 
with  each  pound  of  gas  was  evaluated  by  determining  the  total  equili¬ 
brium  enthalpy  of  the  escaping  gas  (see  Section  2d)  above  the  enthalpy 
of  the  unheated  explosive.  This  quantity  of  heat  was  determined  by 


substituting  the  enthalpy  into  the  following  equation- 

M 

Q  (gas) 

■  w  J  H<T>dM 

whe  rc 

0 

Q(gas) 

=  heat  per  pound  of  escaping  gas,  B/lb 

H(T) 

=  enthalpy  of  gas  above  its  value  at  ambient 

temperature,  B/lb 

dM 

=  increment  of  gas  mass,  lb 

The  quantity  of  heat  released  by  the  burning  of  1  pound  of 
the  high- explosive,  Q^.  is  expressible  as 

Qj.  =  Q(chamber)  +  Q(gas)  +  Q(insulation)  +  Q(misc)  (3) 

Where  Q(misc)  refers  to  the  small  quantities  of  heat  associated 
with  the  heater,  ignitor,  and  kinetic  energy  of  the  escaping  gas. 

In  addition,  a  literature  study  was  conducted  which  revealed 
that  the  quantities  of  heat  released  by  the  slow  decomposition  of  TNT, 
Cyclotol,  75/25  and  RDX  were  900,  675,  and  900  B/lb,  respectively  , 
while  the  quantities  of  heat  released  by  the  explosion  of  TNT,  Cyclotol, 
75/25,  Composition  B,  RDX  and  HMX  were  1944,  2205,  2232,  2304, 
and  2441  B/lb,  respectively  (see  Table  A  1).  The  quantities  of  heat 
released  by  the  explosions  are  considerably  higher  than  those  from 
(lie  slow  decomposition.  Since  the  velocities  of  these  two  types  of 
reactions  bracket  the  burning  velocities  from  the  ARF  experiments, 
it  may  be  expected  that  the  values  of  the  heat  released  by  burning  a 
pound  of  a  given  explosive  would  lie  between  the  values  associated 
with  the  two  reactions.  To  determine  whether  or  not  this  is  true, 
the  heat  released  by  a  pound  of  Composition  B  was  determined  as  a 
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function  of  the  velocity  of  burning.  The  results  are  shown  in  Figure  2. 
The  rate  of  heat  release  from  Composition  B  at  “zero  burning  velocity" 
was  considered  equal  to  that  from  Cyclotol  75/25  during  slow  decom¬ 
position,  It  may  be  noted  that  the  values  are  reasonably  consistent 
and  that  the  heat  released  by  the  burning  Composition  B  reaches  a 
maximum  at  the  relatively  low  velocity  of  about  40  ft/hr. 

The  heat  released  by  burning  PBX  9404  under  atmospheric 
pressure  was  also  measured  and  resulted  in  the  following  measure¬ 
ments:  2223,  2261,  24.63,  2492  and  2736  B/lb.  The  average  of 
these  values  is  2433  B/lb.  Unfortunately  no  values  could  be  found 
in  the  literature  for  the  heat  of  explosion  of  PBX  9404,  Since  the 
major  constituent  of  PBX  94  04  is  HMX,  it  appears  reasonable  that 
the  heat  of  explosion  of  PBX  94  04  would  be  close  to  tnat  of  HMX, 
i.  e. ,  2441  B/lb.  It  may  be  observed  that  the  difference  between 
this  value  and  the  average  value  for  slow  burning  is  small  and  is 
much  less  than  the  experimental  variations.  From  this  it  appears 
that  the  heat  released  by  burning  PBX  9404  is  essentially  indepen¬ 
dent  of  burning  velocity. 

b.  Enthalpy  of  Composition  B  and  HMX 

The  enthalpy  of  a  high-explosive  plays  a  very  significant  role 
in  affecting  its  rate  of  burning.  This  is  due  to  the  fact  that  the  vir¬ 
gin  material  must  first  be  heated  to  temperatures  of  about  500  to 
600*F  before  the  material  will  react  with  any  significant  vigor. 
Unfortunately,  no  enthalpy  data  could  be  found  above  about  300* F. 

Values  of  the  enthalpy  of  Composition  B  and  HMX  may  be  obtained 
from  the  specific  heat  values  given  in  Table  A-l.  Curves  of  enthalpy 
versus  temperature  constructed  from  the  data  from  Composition  B 
and  HMX  are  shown  in  Figure  3,  These  curves  were  extrapolated 
to  the  higher  temperatures  by  assuming  that  the  heats  of  fusion  of 
HMX  and  RDX  are  equivalent  to  that  of  TNT.  Due  to  the  uncertainties 
of  this  extrapolation,  the  extrapolated  portions  of  the  curves  are 
shown  as  dashes.  No  account  was  taken  for  possible  vaporization 
of  the  high-explosive  material.  The  curve  for  Composition  B  was 
taken  from  the  DASA  Report  1234. 

c.  Conductive  Heating  of  Burning  Composition  B 

During  burning,  an  explosive  is  heated  by  radiation  (to  be 
discussed  later)  and  by  conduction  from  the  hot  gases.  At  pressures 
in  excess  of  a  few  hundred  atmospheres,  conduction  is  the  dominant 
source  of  heat  flow  to  the  reacting  explosive.  A  literature  search 
was  conducted  to  find  any  data  which  either  would  indicate  the  con¬ 
ductive  heat  flux  or  could  be  used  to  determine  the  flux.  Fortunately, 
a  substantial  amount  of  strand  burn. data  was  found  in  a  Picatinny 
Arsenal  report  for  Composition  B.  These  data  indicate  the  burning 
velocity  as  a  function  of  pressure  and  were  obtained  by  burning  a 
long  thin  rod  of  the  explosive  in  an  environment  which  was  maintained 
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at  constant  pressure.  This  was  accomplished  by  use  of  a  chamber 
whose  volume  is  large  compared  to  that  of  the  gas  evolved  by  explosion. 
The  velocities  of  burning  of  Composition  B  are  presented  in  Table  I 
for  pressures  ranging  from  136  to  2041  atmospheres.  Since  thermal 
radiation  is  not  a  significant  source  of  heat  flux  at  these  pressures, 
it  is  possible  to  utilize  the  Arrenhius  relationship  and  the  data  gen¬ 
erated  in  the  DASA  Report  1234  to  determine  the  conductive  heat 
flux  to  the  reacting  explosive  at  the  several  pressures.  In  IITRI's 
previous  study,  it  was  found  that  the  thickness  of  the  layer  of  Com¬ 
position  B  undergoing  reaction  was  about  0.  05  inch  or  0.  0042  foot. 

This  layer  consisted  of  foam  whose  density  was  found  to  range  bet¬ 
ween  1/3  to  1/6  of  that  of  the  solid  explosive.  This  was  determined 
by  slowly  heating  a  bath  of  liquid  Composition  B  until  it  buret  into 
flame.  In  the  present  study  the  density  of  the  foam  layer  will  be 
considered  as  l/6  of  that  <vf  the  virgin  Composition  B.  This  yields 
a  foam  density  of  105/6  or  17,  5  lb/ft3.  Tne  Arr.  nhius  relationship, 
which  was  discussed  in  the  DASA  Report  1234,  is  as  follows 

W  =  P^Zexp.  -(E/RgQTf)  (4) 

where 

W  =  rate  of  mass  flow  per  unit  area,  lb/ ft2  hr 
Pf  =  density  of  foam  =  17.5  lb/ft3 

Xf  =  thickness  of  foam  =  0.  0042  ft 

Z  =  frequency  factor  =  1.  14  x  102  2  /hr 

E  =  activation  energy  =  4 7,  500  Cal/mole 

Rro  =  univeral  gas  constant  =  1 .  1 04  Cal/mole  * R 

Tf  »  temperature  of  foam,  *R 

The  values  of  Z  and  E  are  assumed  to  be  independent  of 
temperature  and  equal  to  those  of  Cyclotol  75/25  shown  in  Table  A  1. 
While  the  values  undoubtedly  vary  with  temperature,  the  inconsis¬ 
tencies  of  the  values  in  the  literature  tend  to  overshadow  this  depen¬ 
dence. 


Substituting  the  appropriate  quantities  into  Eq.  4  yields 

W  =  (17.  5)  (0.  0042)  (1.  14  x  102  2  )  exp.  -(47.  500/l.l04Tf)(5) 

=  8.  38  x  1020  exp.  -(43,  025/Tf) 

The  rate  of  mass  flow  may  also  be  determined  from  the  rate 
of  burning  v  by 

W  =  PEV  (6) 

where 

p„  =  density  of  virgin  Composition  B  =  105  lb/ft3 

v  =  velocity  of  burning,  ft/hr 

Eliminating  W  from  Eq.  S  and  6  solving  for  T^  yields 

Tf  =  43,  025/ln(7.  981  x  10lft/V)  (7) 
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TabU  I 

STRAND  BURNING  RATE  OF  COMPOSITION  B 


Pressure. 

Atmospheres 

Burning  R&te, 
ft/hr 

136 

153 

204 

246 

340 

480 

680 

1071 

1361 

2001 

2041 

3126 
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This  indicate*  the  temperature  of  the  foam  as  a  function  of 
the  velocity  of  burning.  The  enthalpy,  H,  of  the  foam  for  each  of  the 
several  velocities  may  be  found  from  the  data  shown  in  Figure  3. 
Following  this  the  rate  of  conductive  heating,  qc>  required  to  main¬ 
tain  each  of  the  several  velocities  of  burning  may  be  expressed  by 
(  H  =  enthalpy  above  that  at  70*  F) 

qc  =  PEH(Tf)V  (8) 

This  represents  the  rate  of  conductive  heating  as  &  function 
of  the  velocity.  The  relationship  between  the  conductive  heating  and 
pressure  may  be  obtained  by  the  velocity-pressure  data  shown  in 
Table  I.  A  plot  of  the  ratio  of  the  conductive  heating  and  the  pressure 
is  shown  in  Figure  4  as  a  function  of  pressure.  A  ratio  was  used  to 
maintain  a  reasonable  scale  which  allows  a  better  appreciation  of  the 
effects  of  pressure.  It  may  be  noted  that  the  conductive  heating  in¬ 
creases  at  a  greater  than  linear  rate  with  pressure  at  least  up  to 
2000  atmospheres  of  pressure.  It  may  also  be  observed  that  at  near- 
atmospheric  pressure  the  value  of  the  conductive  heating  is  about 
25,  000  B/ft2  hr.  This  value  may  be  compared  with  the  value  of 
41,  000  B/ft2  hr  which  was  determined  in  the  previous  study  (Bee 
DASA  Report  1234)  by  the  subtraction  of  the  substantial  radiant  flux 
(experimental)  from  the  total  heat  flux  (theoretical).  Based  on  the 
accuracies  of  the  techniques  used,  the  agreement  between  the  two 
values  is  gratifyingly  good. 

d.  Total  Equilibrium  Enthalpy  of  Gas  Released  by  Burning 

Composition  B,  RDX  and  HMX 

The  temperatures  of  the  gas  from  burning  explosives  are  suf¬ 
ficient  to  cause  appreciable  dissociation  and  ionization.  As  a  result 
of  this  the  specific  heat  increases  significantly  with  temperature  and 
is  not  accurately  treated  as  a  constant  value  over  the  temperatures 
of  concern.  The  total  enthalpy  of  the  gas  evolved  during  the  burning 
of  Composition  B,  RDX  and  HMX  was  computed  as  a  function  of 
temperature  from  900*R  to  6300*R  assuming  equilibrium  conditions. 

This  was  accomplished  by  determining  the  composition  of  the  gas 
and  summing  the  heats  of  formations  and  enthalpies  of  the  various 
species  at  each  of  several  temperatures.  The  composition  of  the 
gas  was  determined  using  the  techniques  illustrated  by  Huff4  and 
Eberhardt.  The  chemical  formulas  of  the  high-explosives  are 
as  follows: 


COMPOSITION  B 

.  605 

moles 

C3H6N6°6  <RDX> 

.  395 

moles 

C?H5N306  (TNT) 

HMX 

1 

mole 

C4H8N8°8 

RDX 

1 

mole 

C3H6N6°6 

1 1 


70.  000 
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It  may  be  observed  that  the  number  of  atoms  in  the  HMX  and 
RDX  molecules  differs  only  by  a  factor.  As  a  result  the  composition 
and  enthalpies  of  the  gases  associated  with  HMX  and  RDX  will  be 
identical  under  conditions  of  equilibrium  and  for  equal  temperatures  and 
equal  pressures.  The  percentages  of  the  gas  constituents  associated 
with  Composition  B  and  HMX  or  RDX  are  summarized  in  Table  II 
and  Table  III  for  a  pressure  of  1  atmosphere.  The  sums  shown  at 
the  bottom  of  each  column  are  not  all  100  percent  due  to  errors  in 
rounding  off  the  numbers.  Calculations  were  also  performed  at  10 
and  100  atmospheres.  These  calculations  predicted  small  shifts  in 
the  composition  of  the  gas.  However,  the  effect  of  these  shifts  on  the 
total  enthalpy  is  quite  small  and  will  be  neglected  in  the  present  study. 


The  total  equilibrium  enthalpy,  H(T),  was  computed  using 
the  equation 

Hm  =  'ir  E  si  [  Hi  +  m 

where  the  factor  of  1.8  is  introduced  to  give  values  of  H(T)  in 
B/lb  rather  than  Cal/gm  and  where 

M  =  effective  molecular  weight  of  gas 

S.  =  moles  of  the  ith  constituent  per  mole  of  gas 

H£  =  enthalpy  (sensible  heat)  per  mole  of  the  ith 

constituent 

(AI-J^.).=  heat  of  formation  of  the  ith  constituent 

The  resulting  values  of  the  total  equilibrium  enthalpies  of 
the  Composition  B  and  HMX  and  RDX  gas  were  referenced  to  a  tem¬ 
perature  of  530"R  (70*F)  and  are  shown  in  Figures  5  and  6  respectively. 
In  both  cases  it  may  be  observed  that  the  total  enthalpy  is  definitely 
not  linear  with  temperature.  Therefore  the  specific  heat,  which  is 
equal  to  the  slope  of  the  curve,  changes  vary  appreciably  with  temper¬ 
ature,  This  is  particularly  true  at  the  high  temperatures.  (The 
specific  heats  used  here  refer  to  the  slopes  of  the  total  equilibrium 
enthalpy  curves  rather  than  the  normal  enthalpy  curves  in  which  the 
gas  constituents  remain  invariant.  ) 


From  the  data  in  Tables  II  and  III,  it  is  a  simple,  exercise 
to  determine  the  gas  constant  and  the  ratio  of  the  "specific  heat"  as 
constant  pressure  to  the  "specific  heat"  at  constant  volume  for  each 
of  the  specified  temperatures.  To  accomplish  this,  the  apparent 
molecular  weight  of  the  gas  was  found  by  taking  a  weighted  average 
of  the  molecular  weights  of  the  several  constituents.  The  individual 
gas  corystant8  may  be  evaluated  by  dividing  the  universal  gas  constant, 
1544 -r i  p  >  by  the  apparent  molecular  weights.  The  ratios  of 
the  ^'specific  heats,  "  k,  may  be  found  from  the  following 

well  known  equation 


k  -  JC  /(JC  -  R  ) 
P  P  g 


(10) 
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The  results  are  summarized  in  Tables  IV  and  V>  where 

J  -  mechanical  equivalent  of  heat  -  ?78  ft=lfc/B 

R  =  gas  constant,  ft-lb/lb  *R 

8 

c.  Specific  Heat  and  Thermal  Conductivity  of  Steel  Chamber 

As  described  in  Section  2  a,  the  steel  chamber  was  used  as  a 
calorimeter  to  measure  the  heat  evolved  by  burning  explosive.  Later 
in  the  report  the  temperature  distribution  along  the  length  of  the  chamber 
will  be  computed  and  compared  to  the  experimental  temperatures. 

This  will  serve  as  one  of  several  checks  on  the  radiation  theory  pre¬ 
sented  in  a  subsequent  section.  Accurate  values  of  the  specific  heat 
of  the  steel  as  a  function  of  temperature  are  required  in  both  cases. 

In  addition  a  knowledge  of  the  thermal  conductivity  is  necessary  for 
the  computation  of  the  temperature  distribution  along  the  length  of  the 
chamber.  Values  of  the  specific  heat  and  thermal  conductivity  of 
ordinary  carbon  steel  found  in  the  literature”  are  listed  in  Table  VI. 

Due  to  the  importance  of  the  specific  heat  for  the  calorimetric  measure¬ 
ments,  three  enthalpy  measurements  of  the  steel  were  made  between 
72°  and  21 2°  F  and  also  between  72*  and  500*F.  These  measurements 
were  made  by  dropping  a  heated  sample  of  the  steel  into  a  water  bath 
contained  in  a  Dewar  calorimeter.  The  temperature  of  the  water 
bath  was  measured  by  a  Beckman  thermometer.  These  enthalpy 

measurements  resulted  in  the  following  average-values  of  specific 

i — .  r 

stl 

C  .  =  .111  B/lb*F,  72  to  21 2* F 

C8'  =  .  116  B/lb*F,  72  to  500*F 

stl 

It  may  be  observed  that  these  average  values  compare  quite 
favorably  with  the  values  shown  in  Table  VI. 

f.  Thermal  Radiation  From  Gas  Evolved  by  Burning  Explosives 

A  significant  portion  of  the  heat  flux  to  a  burning  explosive 
at  low  pressures  (of  the  order  of  tens  of  atmospheres)  is  generally 
due  to  thermal  radiation. 

Thermal  radiation  is  emitted  by  the  individual  molecules  and 
by  the  particles  within  the  gas  stream.  The  gases  evolved  by  TNT- 
based  explosives  such  as  Cyclotols  and  Composition  B  contain  appre¬ 
ciable  quantities  of  Buch  particles.  The  emission  of  radiation  from 
the  gas  molecules  evolved  by  the  burning  of  such  explosives  is  neg¬ 
ligibly  small  compared  to  that  from  the  cloud  of  particles.  This  will 
become  more  evident  later  in  the  report. 

The  first  part  of  this  section  will  be  devoted  to  a  derivation 
of  the  basic  equation  that  prescribes  the  radiant  emission  from  the 
particles  contained  in  a  gas  at  constant  temperature  and  pte  ssure. 
Following  this,  the  equation  will  be  used  to  determine  the  equations 
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Table  IV 


PROPERTIES  OF  GAS  EVOLVED  BY  BURNING  COMPOSITION  B 


Temperature 
”  R 

Molecular 

Weight 

Gas  Constant 
ft  -  lb/lb  °R 

Ratio  of  "Specific 
Heats" 

900 

21.4 

72.  1 

1.  37 

1800 

21.  5 

71.  8 

1.  34 

2700 

21.  5 

71.  8 

1.  29 

3600 

21.  3 

72.  5 

1.  28 

4500 

21.  2 

72.  8 

1.  20 

5400 

20.  6 

75.  0 

1.  11 

6300 

19.  0 

81.3 

1.  06 

Table  V 


PROPERTIES  OF  GAS  EVOLVED  BY  BURNING  RDX  OR  HMX 


Tempi*  r.i  tu  re 
‘  *  R 

Molecular 

Weight 

Gas  Constant 
ft- lb/lb  *R 

Ratio  of  "Specific 
Heats" 

900 

24.  6 

62.  7 

1.  24 

1800 

24.  8 

62.  3 

1.  24 

2700 

24.  7 

62.  5 

1.  23 

3600 

24.  6 

62.  8 

1.  19 

4500 

24.  6 

62.  8 

1.  16 

5400 

22.  6 

68.  3 

1.  08 

6300 

19.  6 

78.  8 

1.  03 
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Table  VI 

n&Ai  rtiiD  THERMAL  CONDUCTIVITY 
OF  PLAIN  CARBON  STEEL 
(DENSITY  =  491  LB/FTj) 


Thermal 

Temperature 

Specific  Heat 

Conductivity 

*F 

B/lb*F 

B/ft*  hr*F 

70 

.  100 

.37.  2 

100 

.  104 

36.7 

200 

.  no 

35.  1 

300 

.116 

■•3.  5 

400 

.  120 

32.  0 

500 

.  124 

30.  3 

600 

.  134 

29.  0 

700 

.  140 

27.  3 

800 

.  159 

26.  0 

900 

.  164 

24.4 

1000 

.  179 

23.  1 

1100 

-  -  -  - 

21.8 

1200 

— 

20.  5 

1300 

-  -  -  - 

19.  0 

Z0 
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describing  the  radiation  emitted  by  a  cylindrical  volume  of  the  gas. 
Thermal  radiation  will  be  divided  into  two  parts. . .  that  arriving 
at  the  circumferential  surface  enclosing  the  cylindrical  volume  of 
gas  and  that  arriving  at  the  two  faces  at  either  end  of  the  gas  volume. 

(1)  Derivation  of  Fundamental  Radiation  Equation 


There  are  four  principal  parts  to  the  radiation  equation. 

These  are: 

(a)  The  emission  of  radiation  from  a  gas  element. 

(b)  The  divergence  of  the  radiation  in  passing  from  an 
element  to  the  target  area. 

(c)  The  attenuation  of  the  radiation  in  passing  from  an 
clement  to  the  target  area. 

(d)  The  orientation  of  the  gas  element  and  target  area 
with  respect  to  each  other. 


Figure  7  represents  a  gas  element  and  an  increment  of  sur¬ 
face  being  irradiated.  The  radiation  emitted  by  the  particles  in  a 
gas  clement  must  be  proportional  to  the  size  and  density  of  the 
particles.  If  equal  masses  of  the  gas  contain  the  same  number 
anH  of  nartirles  the  emitted  radiation  will  ho  proportional 

to  the  gas  density  P  .  In  addition  the  emission  of  radiation  is 
inversely  proportionSl  to  the  wave  length  A  raised  to  a  power  close 
to  1.  For  purposes  of  simplifying  the  expressions,  the  exponent 
will  be  taken  as  1.  Based  on  these  arguments,  the  emission  of 
heat  radiation  is  approximately  equal  to  the  factor  p  p  dR./A  ,  where 
p  is  a  constant  of  proportionality  and  dR  is  the  thiclmess  of  the  gas 
crciTient  in  the  direction  prescribed  by  the  line  between  the  element 


and  the  surface  being  irradiated.  The  amount  of  radiation  emitted 
within  a  given  increment  of  wave  length  is  proportional  to  the  product 
of  monochromatic  a  teal  radiant  intensity,  N^.  °  and  the  increment 
of  wavelength.  dA  As  a  result  of  geometrical  considerations,  only 
a  fraction  of  the  omitted  thermal  radiation  is  incident  on  the  surface 
being  radiated.  This  fraction  is  proportional  to  the  product  of  the 
component  of  the  vector  representing  the  area  of  the  gas  element  in 
view  of  the  irradiated  surface.  dA  cos  0,  and  the  solid  angle  sub¬ 
tended  by  the  surface  being  irradiated,  dA  cos$/R2  .  A  portion 
of  this  radiation  is  intercepted  by  the  cloud  of  the  intervening  particles. 
To  determine  the  amount  that  is  intercepted,  consider  a  collimated 
beam  of  intensity  I  and  of  wave  length  A  passing  through  an  element 
having  a  thickness  of  dR.  The  intensity  of  this  beam  will  decrease  in 
the  following  fashion 


dl  =  -  V5  p  dR/A 
g  g 


(ID 


Since  the  pressure  and  temperature  of  the  gas  are  considered 
uniform  P  as  well  as  ^  and  A  are  constants.  Integration  of  this 
exprescion  with  I  =  I*  at  R  =  0  gives 

l/V  =  exp.  -(3  p  R/A) 
g  g 


(12) 
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This  indicates  that  the  fraction  of  the  radiation  which  passes 
through  a  gas  element  of  thickness  R  is  exp.  -(£  P  RAJ*  As  a  result 
the  radiation  which  is  emitted  from  a  gas  elemem  with  wave  lengths 
between  A  and  A  +  dA  and  is  incident  on  the  surface  may  be 
described  as  follows 


dq  *  P  - 
r  g 


Nn,(T  )dAdAo  cosB-^-^Ll  exp-{P  p  R/A) 
is  a  a  a  rs  2  K  K 


If  dV  represents  the  volume  of  the  gas  element,  then  dV  =  dA  dz, 
and  since  dz  =  dR  cos  0,  8  8 

dA  =  dV  /dR  cos  9  (14) 

g  g 

Substitution  of  this  into  Cq.  13  yields 


-  .^SPgN%.:i!Ig)  exp  -(P  p  R/A)dAdV  dA 


where 

q  »  radiation  incident  on  target  area,  B/ft2  hr 

fig  ■  constant  of  proportionality 

Pg  =  density  of  gas.  lb/fts 

A*’  =  wave  length,  cm 


A5  [exp  (C/Aiy-l]' 


B/ft2  hr  cm4 


absolute  temperature  of  gas  or  particles,  *R  „  4 

first  radiation  coefficient  -  1.86?  x  10"°  - 

second  radiation  coefficient  -  2.58  cm  *R  ft2  hr 

slant  distance  between  gas  element  and  target  surface,  ft 
volume  gas  olement,  ft* 
incremental  area  of  target  surface,  ft2 


It  may  be  noted  that  the  temperature  of  the  gas  and  particles  are 
used  interchangeably.  This  is  a  result  of  an  analysis  of  the  heat  trans¬ 
fer  to  the  particles  by  Schack'  which  indicates  that  the  temperatures 
of  the  gas  and  particles  are  essentially  equal,  due  to  the  small  dimen¬ 
sions  of  the  particles. 

To  simplify  the  expression  for  dqf,  AT  will  be  replaced  by  Aq. 
Also  the  ideal  gas  law  will  be  considered  valid,  i.  e. , 

P 


P  =  gas  pressure,  lb/ft2 

R|  =  gas  constant,  ft- lb/lb  *R 
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Integrating  the  part  of  Eq.  15  containing  A  and  substituting 
the  relationships  of  Eqs.  16  and  l?  gives 


P„P0 
8  g 


R  exp-(PgPgRA)  dX  = 


>2 C.exp  -(P  P  R/R  A  ) 

J  ^  8  ~8 _ a_°- 


P  P  RdX 

_JL_2 _ ? 


[exp[C2/*o)-l] 


08) 


=  T4  F 
g 


P  P  R 

-V- 


g 


The  function  F  is  presented  in  Figure  8  with  pg^Rg  represented 
Substitution  of  the  above  relationships  into  Eq.  ft  yields 


by  Pg.  Substitution  of  the  above  relationships 

dq  =  T4  -S2L*  F(P'P  R)dVdA 
g  g  g  8 


09) 


This  is  the  basic  equation  describing  the  radiant  flux  that  is  emitted 
by  an  element  of  gas  at  uniform  temperature  and  pressure  and  is  inci¬ 
dent  on  an  area  dA.  The  next  section  will  apply  this  equation  to  the 
specific  problems  of  computing  the  radiant  fluxes  from  a  cylindrical 
volume  of  gas. 

(2)  Radiation  Passing  Through  the  Circumferential  Boundary  of 
a  Cylindrical  Volume  of  Gas 

The  thermal  radiation  passing  through  the  boundary  at  the  sides 
of  a  cylindrical  volume  of  gas  generally  represents  a  significant  heat 
loss  from  the  gas.  A  knowledge  of  this  loss  is  of  key  importance  in 
the  analysis  of  the  behavior  of  burning  explosives,  which  will  be 
covered  in  a  later  section.  This  section  will  be  devoted  to  a  derivation 
of  an  equation,  that  will  be  used  later  in  the  report,  to  determine  the 
effects  of  this  heat  loss. 


A  diagram  of  a  volume  and  surface  element  of  a  cylindrical 
volume  of  gas  is  shown  in  Figure  9.  In  the  derivation,  the  temperature 
and  pressure  of  the  gas  will  be  considered  as  constant  throughout  its 
volume. 
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From  trigonometric  relationships 

»/•> 

R  =  (a2  -  2ar  cost  +  r2  +  z'2  )  ' 

cos$  =  (a  -  r  cost  )/R 

AV  a  rAtArAz 

8 

AA  =  a  At 1  Az ' 


(20) 

(21) 

(22) 

(23) 


Substitution  of  these  relationships  into  Eq.  19  yields 

Aq'  =  T*  -fe-.-. r-5..08J  ]  F((3  1  P  R)  rAt  ArAzaAt  'Az' 
8  R«  g  8 


(24) 


The  radiation  which  is  emitted  by  a  disc-shaped  element  of  the 
gas  and  which  passes  through  the  circumferential  boundary  of  the  gas 
volume  may  be  determined  by  integrating  Eq.  24  as  follows 


Aq;  =  4«T*  i 


I  z  a  tt 

\  f  C  f  *■  -■  ^  F(P  '  P  R)  ar  dt  dr  dz' 

J  J  J  R4  8  8 


o  o  o 


(25) 


£  -z  a  tr 

o  o  o 


R4 


-\ 

F(P  'P  R|  a rd»  dr  dr' 

8  8  J 


In  order  to  simplify  this  expression,  let  the  function  P  repre¬ 
sent  the  following 


a  w 


P,  (z,  a,  (3  '  P  )  =  4"a 
1  g  8 


o  o  o 


<a  ~  r..?°P-4.>  F(P'P  R)  rd«  dr  dz' (26) 
R4  8  8 
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Then  Eq.  25  becomes 

AqJ.  =  T*g  |p(a.a,PjPg)  +  P(*  -z,  a,  P^Pg)  J  Az  (27) 

This  represents  the  rate  of  heat  loss  from  a  gas  disc  which  passes 
through  the  circumferential  boundary  of  the  gas  volume  (see  Figure  9). 
Several  plots  of  Pj(z,  a=  \>  6ft,  PgPg)  are  shown  as  functions  of  *  in 
Figure  10  for  several  values  of  PgPg.  It  may  be  observed  that  as  P|F_ 
increases,  P[  approaches  its  maximum  value  at  decreasing  values  of  nt. 

In  the  limit,  where  the  gas  becomes  completely  opaque,  the  rate  of  heat 
loss  would  tend  to  approach  the  radiation  emitted  from  the  edge  of  the 
gas  disc.  The  area  of  this  edge  is  2iraAz,  From  the  Stefan- Boltzmann 
Equation 

limAq'  —>■  2iraAzaT4  as  P'P  — «  (28) 

r  8  g  g 

where  o'  =  Stefan- Boltzmann  Constant  =  .173  x  10*°B  ft*hr*R4, 

Comparison  of  this  expression  with  the  right  side  of  Eq.  27  shows 
that  P]  approaches  tra<r  or  (3.  1416)  (0,  1667)  (0.  173  x  10"®)  =  0.  906 

x  10" 9.  It  may  be  observed  that  this  value  is  approached  quite  closely 
for  PgPg  =  30.  This  completes  the  derivation  of  the  equation  pre¬ 

scribing  the  radiant  heat  loss  from  a  disc  shaped  element  of  gas  through 
the  circumferential  boundary  enclosing  a  gas  volume  of  uniform  temperature. 

(3)  Radiation  Passing  Through  the  Faces  of  a  Cylindrical  Gas  Volume 

The  thermal  radiation  passing  through  the  boundaries  at  either 
end  of  a  gas  volume  represents  an  important  source  of  heat  loss  by 
the  gas.  It  also  represents  a  significant  source  of  heating  of  the  re¬ 
acting  explosive  since  the  explosive  was  situated  at  one  end  of  the  gas 
volume  in  experiments  to  be  discussed  later  in  this  report.  This 
section  will  be  devoted  to  a  derivation  of  an  equation  that  will  prescribe 
this  radiant  flow  for  various  gas  pressures  and  temperatures. 

A  diagram  of  a  volume  element  and  a  surface  element  of  a  cylin¬ 
drical  volume  of  gas  is  shown  in  Figure  11.  It  should  be  observed  that 
each  gas  clement  is  encompassed  by  circles  drawn  from  a  point  lying 
directly  above  the  surface  element.  As  was  true  on  the  previous 
analysis,  the  pressure  and  temperature  will  be  considered  as  constant 
throughout  the  gas.  From  trigonometric  relationships 


cosb 

=  z/R  . 

-  (r*  .  .«  )‘/2 

(29) 

R 

(30) 

AV 

8 

AA 

=  6  (r',  r)  rArAz 

(31) 

-  r'A«  'Ar' 

(32) 

( 

=  2  arc  cos  (r*  +  r2  -  a*  }/  2rr' 

(33) 

28 
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Fig.  11  RADIATION  RECEIVED  BY  END  FACE 
OF  CAS  VOLUME 
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Substitution  of  these  relationships  into  Eq.  19  gives 

Aq"  =  T4  F{P'P  R)  |(r',  rJrArAzr'Ac'Ar' 

r  8  r*  g  g 

From  Figure  11,  it  may  be  noted  that 

4  =  Zff  for  r  £  a-r* 

o<{<2tr  for  a-r*  <  r  <■  a+r' 

4  =  o  for  r  a+r 1 


(34) 


a-r' 


Aq"  =  T* 

r  g 


fi 


F((3  'PR)  2tt  r  dr 
g  g 


(35) 


a+r' 


F(p  '  P  R)  4  (r\ r)  r  dr  >  2irr'dr'Az 
g  g 


A'  R* 

This  equation  will  be  simplified  by  defining  a  function  Sj  as 


S  (z.  a,  P  ’  P  ) 
1  g  g 


-5—  F(P  1  P  R)  2tt  r  dr 
R4  «  8 


(36) 


a  tr' 


r — 

J  ,  R4 


F(P  '  P  R)4  (r r)  rdr  f  ’  2tr  r'dr' 
g  g 


a-r' 

Then  Eq,  35  becomes 


Aq"  =  T4  S.(z,a,  P'P  )Az 
1  r  g  1  g  g 


(37) 
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This  represents  the  rate  of  heat  loss  from  the  element  of  gas  to 
the  end  surface  of  the  gas  volume  which  is  at  a  distance  z.  The  loss 
to  the  other  end  of  the  gas  volume  is 


Aq”  =  T*g  Sj  (1  -z,a,pMPg)Az  (38) 

The  total  radiant  flux  to  the  end  surfaces  may  be  found  by  inte¬ 
grating  Eq.  37  and  38  with  respect  to  z.  Several  plots  of  Sj  as  functions 
of  z  are  shown  in  Figure  12  for  several  values  of  PgFg>  It  may  be 
noted  that  Sj  increases  with  PiF_.  This  should  happen  since  the  gas 
becomes  increasingly  opaque  witrn  larger  values  of  f3iPg,  In  addition 
it  may  be  expected  from  the  Stefan-Boltzmann  equation  that 


lim  q» 


lim  T4 
g 


ita2o-T4  as  P'P  —y 00 
8  g  8 


(39) 


This  relationship  was  found  to  be  valid  by  integration  of  the 
function  over  sufficiently  large  values  of  z  and  (3iP„.  This  completes 
the  presentation  of  the  basic  data  and  formulas.  77)6  subsequent 
sections  will  be  devoted  to  a  discussion  of  experimental  studies  of 
burning  high-explosive  and  an  application  of  these  data  and  formula 
to  explore  the  effects  and  behavior  of  burning  explosives. 
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i.  BURNING  OF  COtoPOSiiION  b  PBX  9401  UNDER 

PRESSURES  OF  1  TO  3  ATMOSPHERES 

Composition  B  and  PBX  9404  burn  at  speeds  between  about  2  and 
20  ft/hr  under  pressures  from  1  to  3  atmospheres.  Due  to  these  low 
rates  of  burning,  it  is  possible  to  make  a  number  of  reliable  measure¬ 
ments  that  are  not  possible  at  the  higher  burning  rates.  Furthermore, 
a  knowledge  of  the  behavior  of  burning  explosives  under  low  pressures 
is  necessary  to  predict  their  behavior  after  being  ignited  in  environ¬ 
ments  of  near-atmospheric  pressures.  A  number  of  tests  were  con¬ 
ducted  in  this  program  to  determine  the  behavior  of  burning  explosives 
under  small,  slowly- increasing  pressures,  In  these  experiments,  the 
temperature  and  irradiation  were  measured  in  the  gas  stream.  Additional 
experiments  were  conducted  to  determine  the  heat  released  by  the 
burning  (see  Section  2a  ).  The  first  part  of  this  section  will  describe 
the  experiments  and  summarize  the  resulting  data.  The  second  part 
of  this  section  will  develop  an  analytical  description  of  the  burning  and 
its  effects,  and  will  compare  predictions  with  experimental  data, 

a.  Experimental  Studies  of  Burning  Explosives  Under  Low  Pressure 

PBX  9404  burns  approximately  twice  as  fast  as  Composition  B 
and  results  in  considerably  higher  flame  temperatures.  The  high 
temperatures  of  the  gas  stream  are  generally  beyond  the  temperature 
range  of  ordinary  thermocouples  except  at  distances  in  excess  of 
2  feet  from  the  PBX  9404  and  under  conditions  of  slow  burning,  i.  e. , 
pressures  of  about  1  atmosphere.  As  a  result  of  the  high  temperatures, 
the  experimental  studies  of  PBX  9404  were  limited  to  the  unconfined 
burning  in  the  "open"  chamber  indicated  in  Figure  1. 

The  experimental  studies  of  the  burning  of  Composition  B  were 
conducted  under  conditions  in  which  the  evolved  gas  was  "confined" 
as  well  as  "unconfined".  The  confinement  was  adjusted  to  vary  the 
rate  of  increase  of  the  gas  pressure  and  of  the  burning  velocity.  The 
first  part  of  this  section  will  discuss  the  experimental  studies  of  the 
"confined"  and  "unconfined"  burning  of  Composition  B.  This  will  be 
followed  by  a  discussion  of  the  experimental  studies  of  the  "uncon- 
fined"  burning  of  PBX  9404. 

b.  Confined  and  Unconfined  Burning  of  Composition  B 

These  studies  are  a  continuation  of  those  presented  in  the  DASA 
Report  1234.  Except  for  minor  differences,  the  tubular  chamber  in 
which  most  of  the  high  explosive  was  burned  (see  Figure  13),  was 
essentially  the  same  as  that  used  previously  (see  DASA  Report  1234). 

A  ceramic  liner  was  used  to  protect  the  thermocouple  lead  wires,  to 
provide  a  means  of  mounting  an  electric  heater  and  heat  meters,  and 
to  deter  the  conduction  of  heat  to  the  sides  of  the  cylinder  of  explosive. 

To  prevent  the  escape  of  gas,  an  aluminum  gasket  was  used  between 
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the  face  at  each  end  of  the  steel  tube  and  the  cover  plate.  Two  rubber 
gaskets  were  used  to  seal  the  hole  through  the  bottom  cover  plate.  This 
hole  was  provided  for  the  several  thermocouple  and  electric  wires.  A 
small  hole  was  drilled  through  the  top  cover  plate  to  vent  the  gas  and 
to  allow  a  slow  increase  in  the  pressure  and  burning  rate.  The 
Composition  B  was  heated  with  a  122  watt  heater  (see  Figure  14)  for 
15  minutes  prior  to  ignition.  Ignition  was  accomplished  by  passing 
10  amperes  through  a  0.  56  ohm  resistance  wire  embedded  across 
the  face  of  the  cylinder  of  Composition  B 

The  thermal  irradiation  at  various  locations  in  the  gas  stream 
was  measured  by  the  heat  meter  shown  in  Figure  14.  The  sensing 
element  consisted  of  a  0.  02- inch  thick  disk  of  Nichrome.  The  disk 
was  kept  as  thin  as  possible  to  ensure  a  more  rapid  response  time. 

Both  surfaces  of  the  disk  and  the  upper  face  cf  the  aluminum  cylinder 
were  blackened  with  stoveblack  so  that  their  values  of  emissivity 
and  absorptivity  approached  1  and  were  nearly  equal.  Under  this 
condition,  it  is  possible  to  determine  the  radiation  incident  on  the 
upper  face  of  the  disk  from  the  following  relationship  (see  previous 
DASA  Reports) 


«r 


=  20‘T4 


m 


(40) 


where 

q"  =  irradiation  in  B/'hr  ft* 

T  n 

o-  -  =  Stefan- Iloltzman  Constant  -  0.  173  x  10”  B/hr  ft2  *R4 

Tm  =  absolute  temperature  of  disk  in  *R 

Gas  temperatures  were  measured  both  by  bare  and  shielded 
28-gage  Chromel- Alumel  thermocouple  wires.  No  apparent  differences 
could  be  detected  between  the  data  collected  with  the  two  types  of 
thermocouples.  However,  it  is  not  unlikely  that  the  measured  tempera¬ 
tures  are  slightly  low  since  shielding  is  only  partially  satisfactory 
under  transient  conditions.  Also,  the  gas  flow  was  not  sufficient 
to  produce  high  convective  heat  transfer  to  the  junction  relation  to 
the  radiant  losses  from  the  thermocouples  to  the  relatively  cool  walls 
of  the  assembly.  An  analysis  indicated  that  the  readings  could  be  low 
as  much  as  200* F. 

Gas  pressures  were  measured  by  a  catenary  type  pressure 
gage  over  the  range  1  to  27  atmospheres.  Due  to  the  excessive  heat, 
the  gage  was  mounted  in  a  15-cubic-inch  pressure  vessel  which  was 
connected  to  the  chamber  by  means  of  a  6-inch  length  of  l/4-inch 
copper  tube.  The  rate  at  which  the  explosive  burned  was  determined 
by  locating  thermocouples  along  the  axis  of  the  cylinder  of  explosive. 

A  rapid  increase  in  the  temperature  of  each  thermocouple  indicated 
the  passage  of  the  burning  front. 
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Preliminary  heating  of  the  surface  of  the  Composition  B  was 
required  to  produce  a  reasonably  vigorous  burning  corresponding  to 
velocities  of  2.  5  ft/hr  or  higher.  Very  sluggish  burning  occurred 
in  tests  in  which  there  was  insufficient  heating  of  the  surface  of  the 
Composition  B  prior  to  ignition.  The  reason  for  this  will  become 
apparent  during  the  discussion  of  the  theoretical  analysis  of  the 
factors  affecting  the  burning.  In  each  test*  the  Composition  B  com¬ 
menced  to  burn  with  a  velocity  of  about  2.5  ft/hr.  The  velocity  re¬ 
mained  at  about  this  level  provided  there  was  sufficient  venting  of 
the  gas.  that  is.  when  the  ratio  of  the  cross-sectional  area  of  the 
vent  hole  to  that  of  the  surface  of  the  burning  explosive  was  in  excess 
of  about  0.  006.  With  an  area  ratio  between  0.  003  and  0.  006,  the 
velocity  increased  slowly  for  a  period  up  to  6  minutes.  The  tempera¬ 
ture  and  l.,e  pressure  of  the  gas  rose  gradually  during  this  period. 

When  a  pressure  of  about  2  to  3  atmospheres  was  reached,  the  rate 
of  burning  and  pressure  increased  very  rapidly  and  either  lifted  or 
ejected  the  top  cover  plate.  Curves  of  typical  temperatures  and 
pressures  of  the  gas  are  shown  in  the  DASA  Report  1234. 

A  summary  of  the  temperatures  of  the  gas  at  various  distances 
from  the  surface  of  the  burning  Composition  B  is  shown  in  Figures 
15  to  19  for  several  rates  of  burning.  It  may  be  observed  that  the 
temperature  decreases  quite  rapidly  with  distance  from  the  burning 
surface  of  the  explosive.  In  addition  it  should  be  observed  that  the 
temperatures  are  greater  with  the  higher  rates  of  burning. 

The  temperatures  of  the  heat  meters  are  shown  in  Figures 
20  to  23  as  functions  of  distance  from  the  burning  surface  for  sev¬ 
eral  rates  of  burning.  It  may  be  noted  that  the  meter  temperature 
increases  according  to  the  proximity  of  the  meter  to  the  explosive 
and  also  with  the  rate  of  burning.  This,  of  course,  ;  what  one  would 
expect  and  is  consistent  with  the  temperature  distributions  of  the  gas. 

As  was  mentioned  earlier,  the  size  of  the  vent  hole  was  varied 
to  determine  its  effect  on  the  pressure  buildup  within  the  chamber,  and 
to  acquire  data  at  higher  burning  rates.  With  sufficiently  large  vents, 
the  Composition  B  burned  to  completion  with  only  a  small  buildup  of 
pressure,  i.  c. ,  less  than  1  atmosphere.  On  the  other  hand,  exces¬ 
sively  small  vents  resulted  in  very  appreciable  pressure  buildups 
which  often  broke  or  permanently  lengthened  the  rods  holding  the  cover 
plates,  and  thereby  caused  ejection  of  the  top  cover  plate.  The  re¬ 
mainder  of  the  chamber  always  remained  intact.  The  effect  of  the 
size  of  vent  may  be  measured  according  to  whether  it  was  or  was  not 
sufficient  to  prevent  such  a  pressure  buildup.  A  summary  of  the  dia¬ 
meters  of  the  vent  hole  which  did  and  which  did  not  cause  an  appreciable 
pressure  buildup  (i.  e. ,  of  about  1  atmosphere)  are  shown  in  Figure  24. 
It  may  be  noted  that  there  is  no  sharp  demarcation  between  the  dia¬ 
meters  of  the  vent  hole  for  the  two  conditions.  Thus  while  diameters 
as  small  as  0.  222  inches  resulted  in  uneventful  burning,  holes  as  large 
as  0.  304  inches  resulted  in  a  noticeable  increase  of  pressure  and 
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Support  wires  imbedded 


ELECTRIC  HEATER 


Full  Scale 


Surface  blackened  with 


HEAT  FLOW  METER  Full  Scale 


Fig.  14  ELECTRIC  HEATER  AND  HEAT  FLOW  METER 
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Fig.  15  GAS  TEMPERATURE  (BURNING  RATE  COMPOSITION  B,  2.  5  FT/HR) 
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Tests  in  which  there  was 
no  abrupt  pressure  buildup 
and  pressures  remained 
between  l-2atmos. 


Tests  in  which  there  was 
an  abrupt  pressure  rise 
causing  ejection  of  top 
cover  plate  (see  Fig.  13) 


Fig.  24  EFFECT  OF  VENT  HOLE  ON  BURNING  OF 
COMPOSITION  B 
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caused  ejection  of  the  top  cover  plate.  Therefore!  the  critical  dia¬ 
meter  of  the  vent  hole  associated  with  the  burning  of  Composition  B 
in  the  chamber  shown  in  Figure  13  lies  somewhere  between  0.  222 
and  0.  304  inches.  This  topic,  as  well  as  previous  ones  in  this  section, 
will  be  returned  to  later  in  this  section. 

A  number  of  chamber  temperatures  were  acquired  during  the 
tests  in  which  Composition  B  was  burned  in  an  "unconfined"  manner 
in  the  chamber  shown  in  Figure  1.  The  Composition  B  cylinders  were 
4  inches  in  diameter  and  4  inches  in  length.  Temperatures  were 
measured  along  the  length  of  the  chamber  and  serve  as  one  of  several 
means  of  checking  the  theoretical  work.  A  temperature  profile  along 
the  length  of  the  chamber  is  shown  in  Figure  25,  for  the  case  of 
"unconfined"  burning  of  Composition  B.  These  data  were  obtained 
with  7  thermocouples  located  along  the  length  of  the  chamber  by 
averaging  the  temperatures  from  6  tests.  The  burning  rate  was  2.  5 
ft/hr  and  remained  very  uniform  during  each  of  the  testB.  It  may  be 
observed  that  the  chamber  reached  its  maximum  temperature  at 
about  0.5  feet  from  its  base.  This  was  due  to  the  fact  that  the  lower 
portion  of  the  chamber  was  thermally  protected  at  the  early  times  by 
the  presence  of  the  cylinder  of  Composition  B  and  at  the  later  times 
by  the  layer  of  Adiprene  used  to  bond  the  Composition  B  to  the  chamber. 

c.  Unconfined  Burning  of  PBX  9404 

Cylinders  of  PBX  9404  of  4  inch  diameter  and  5  inches  long  were 
also  burned  in  the  chamber  shown  in  Figure  1  to  determine  the  heat 
evolved  by  the  burning.  These  results  are  given  in  Section  2  a.  The 
teat  procedure  was  identical  to  that  for  the  "confined"  burning  of 
Composition  B.  During  the  course  of  this  work,  a  number  of  chamber 
temperatures  were  procured.  Unfortunately,  due  to  the  very  high 
temperatures  of  the  gas,  measurements  of  the  gas  temperature  could 
only  be  made  near  the  exhaust  end  of  the  2  foot  chamber.  To  over¬ 
come  this  problem,  the  chamber  was  lengthened  to  4  feet.  Three 
tests  were  conducted  with  the  longer  chamber.  In  each  of  these  tests, 
as  well  as  in  the  earlier  tests,  the  PBX  94  04  burned  very  uniformly 
at  a  rate  of  5  ft/hr.  Gas  and  chamber  temperatures  from  the  4-foot 
chamber  are  presented  in  Figure  26.  At  a  distance  of  2  feet  from  the 
burning  surface,  the  temperature  of  the  gas  from  the  unconfined 
burning  of  Composition  B  (see  Figure  15)  was  only  about  960* F  as 
compared  to  about  252Q*F  for  PBX  9404.  This  is  due  to  the  faster 
burning  rate  of  PBX  9404,  Also  PBX  9404  reacts  almost  completely 
at  these  low  burning  rates  while  Composition  B  does  not. 

Summaries  of  the  temperatures  along  the  lengths  of  the  2-foot 
and  4-foot  chambers  are  presented  in  Figures  27  and  28.  It  may  be 
observed  that  the  shapes  of  the  curves  are  similar  to  that  for  Com¬ 
position  B.  The  fact  that  the  maximum  temperatures  were  not  at  the 
base  of  the  chamber  is  due  to  the  protection  afforded  the  lower  5  inches 
of  the  chamber  wall  by  the  explosive  and  by  the  0.  050  inch  layer  of 
Adiprene  cement  used  to  bond  the  high- explosive  to  the  chamber. 
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Fig.  25  CHAMBER  TEMPERATURES  (BURNING  RATE  OF  COMPOSITION  B.  2.  5  FT/HR) 
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Fig.  27  CHAMBER  TEMPERATURES  (B  JRNING  RATE  PBX  9404,5.  0  FT/HR) 
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d.  Theoretical  Analysis  of  the  Effects  of  Burning  High- Explosives 

Under  Low  Pressure 

The  principal  parameters  affecting  the  burning  rate  of  high- 
explosives  such  as  Composition  B  and  PBX  9404  are  the  temperature 
and  the  amount  of  the  explosive  undergoing  reaction.  An  indication 
of  the  temperature,  thickness  and  density  of  the  reacting  layer  was 
gained  from  the  Composition  B  experiments  noted  earlier.  These 
measurements  indicated  that  the  layer  is  in  the  form  of  a  foam  which 
has  a  thickness  of  approximately  0.  05  inches.  The  temperature  of  the 
foam  during  the  burning  of  Composition  B  under  near  atmospheric 
pressure  was  about  540*F.  Throughout  this  study  it  will  be  assumed 
that  the  amount  of  reacting  material  (or  the  thickness  and  density  of 
the  foam)  is  constant.  It  is  obvious  that  the  temperature  of  the  foam 
must  be  maintained  or  increased  to  perpetuate  the  burning.  This 
can  only  be  brought  about  by  a  net  heat  flux  to  the  foam  layer.  As  a 
result,  this  study  will  be  initiated  by  a  consideration  of  the  effect  of 
the  heating  this  layer  on  the  burning  rate  of  the  explosive. 

As  was  mentioned  earlier,  the  principal  sources  of  heat  flow 
into  the  layer  of  foam  are  conduction  and  radiation  from  the  hot  gas 
stream.  The  rate  at  which  heat  is  carried  away  from  this  layer  is 
approximately  equal  to  the  product  of  the  enthalpy  of  the  explosive 
in  the  layer  (above  that  of  the  material  entering  the  layer)  times  the 
r~tc  -t  which.  is  svolvsd.  In  this  ^r^lysiSi  the  hc»t  lest  by 

convection  and  conduction  into  the  solid  explosive  beneath  the  foam 
is  considered  to  be  reabsorbed  by  the  moving  layer.  The  results 
of  this  discussion  may  be  expressed  in  the  following  mathematical 

dH(T  ) 

Vi  — ar —  +  Vf<H(Tt'  •  H(V>Z  •”'p-<E/Vrf  > 1411 


rate  of  heat  flow  into  layer  due  to  conduction  and 
radiation,  B/ ft*  hr 

time,  hr 

density  of  foam  =  17.  5  lb/ ft 
thickness  of  layer,  =  0.0042  ft 
enthalpy  of  foam,  B/lb 

enthalpy  of  high-explosive  entering  foam  layer 
(incoming  explosive),  B/lb 

frequency  factor,  l/hr 

activation  energy  of  explosive,  cal/mole 
universal  gas  constant,  cal/mole  *R 
temperature  of  foam  material,  *R 


formula 

qx(t)  = 

where 

qT(t)  = 


Pc 

* 

xf 

H(Tf)  = 
H(Tfo)= 


Z 

E 
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Throughout  this  report  all  values  of  enthalpy  are  referred  to  a 
temperature  of  70* F. 

For  steady- state  conditions,  i.  e. ,  steady  uniform  burning. 


dH(Tf) 

-ar1* =  0  (42) 

Therefore  Eq.  41  becomes 

qT  =  Pf*f<H<Tf)  -  H(Tfo))Z  exp  “(E/RgQT£)  (43) 

The  velocity  of  burning,  V,  may  be  expressed  as 

PEV(  Tf  )  =  p^Z  exp  -(E/R  Tf)  (44) 


where 

PE  =  density  of  original  explosive  lb/ ft3 
V(  T£  )=  velocity  of  burning,  ft/hr 

Substitution  of  this  expression  into  Eq.  43  yields  the  following 
relationship  for  steady  burning 

qT  =  PEV<  Tf  )(H(Tf)  -  H(Tfo))  (45) 


This  is  the  same  relationship  expressed  in  DASA  Report  1234. 


The  dominant  parameter  in  this  equation  is  the  temperature  of  the 
foam,  Tf.  As  may  be  noted  in  Eq.  44,  small  changes  of  this  tempurature 
will  produce  very  large  changes  in  the  burning  velocity  in  comparison 
to  the  changes  in  the  enthalpy  of  the  foam.  The  result  is  that  for 
steady-state  burning,  the  velocity  of  burning  is  approximately  propor¬ 
tional  to  the  heat  flux  entering  the  explosive. 


The  next  topic  of  concern  is  the  enthalpy  of  the  gas  just  after 
reaction  of  the  vapors  released  by  the  foam.  The  enthalpy  of  the  vapors 
prior  to  reaction  is  approximately  equal  to  H(Tf),  (see  Figure  3).  The 
burning  of  a  pound  of  the  vapor  releases  a  quantity  of  heat  Qg  (see 
Figure  2)  which  is  added  to  the  heat  or  enthalpy  already  in  the  vapor. 

By  conduction  a  certain  amount  of  this  heat  is  returned  to  the  foam 
within  a  few  mean  free  paths  of  the  molecules.  Since  this  heat  loss 
occurs  so  close  to  the  foam,  the  enthalpy  of  the  gas  shortly  after 
leaving  the  foam  is 


H<V 


H(Tf)  +  Qe 


q 


(46) 
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where 

H<Tg) 

H(T{) 

°E 

% 

P' 

g 

W 


enthalpy  of  gas.  B/lb 

enthalpy  of  foam  at  temperature  Tf,  B/lb 

heat  released  by  burning  of  1  pound  of  explosive,  B/lb 

rate  of  conductive  heating  per  atmosphere  of  pressure, 
B/ft2  hr  atmos. 

gas  pressure,  atmos. 

rate  of  gas  flow,  lb/ ft2  hr 


The  terms  H(T^)  and  are  the  dominant  terms  in  this  equation 
and  both  will  tend  to  increase  with  the  burning  velocity.  As  a  result 
of  this,  the  gas  temperature  will  increase  with  the  velocity  of  burning. 
The  increase  of  gas  temperature  may  be  as  much  as  200  to  300* F. 


Since  the  gas  temperature  may  be  as  high  as  6000* R,  the  radiant 
heat  loss  from  the  gas  is  appreciable.  The  rate  of  radiant  heat  loss 
may  be  determined  by  the  equations  derived  in  Section  2..  Let  us 
consider  a  gas  element  of  radius  a  and  thickness  Az  that  is  located 
at  a  distance  of  z  from  one  end  of  a  cylindrical  volume  of  gas  having 
a  temperature  T(z),  and  a  length  £.  From  Eq.  27,  37  and  38,  the  rate 
of  radiant  heat  loss  from  this  element  is 


„  .  ...  P'P)  +  P,  (£-z,  a,  P'P  )  +  S  (z,  a,  p'F  ) 

r  g  L  1  gg  1  gg1  gg 


=  T<j3)  j^fz,  a.  P’P 


(47) 


where  P  and  T  (z)  are  the  pressure  and  temperature  of  the  gas. 
8  g 


This  relationship  presumes  that  the  radiant  flux  between  the 
several  gas  elements  cancel  out.  This  is  considered  to  be  a  good 
first  order  approximation.  The  functions  Fj  and  Sj  are  defined  by 
Eq.  26  and  36  and  are  associated  with  the  heat  losses  through  the 
circumferential  surface  and  the  two  ends  enclosing  the  cylindrical 
gas  volume.  The  expression  within  the  brackets  will  be  simplified  by 
setting  it  equal  to  Uj  (z,  a,  1,  P'  P  ).  The  resuit  of  this  substitution  is 

g  g 

Aqr  =  T4g(z)U1(z,a,  £,p^Pg)Ae  (48) 

Since  Aqr  is  equivalent  to  -ffa2  Azpg(Tg)— ' for  a  particular 

increment  of  gas,  the  rate  of  temperature  change  of  this  gas  element 
as  it  moves  away  from  the  explosive  is 


*a2  Azp  c  (T  )-^p. 
g  P  g  d* 


"  T'gU  j  (2>  a,  *,P-Pg)  Az 


(49) 
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The  total  derivative  of  the  temperature  may  be  written  as 


(50) 


Since  the  rate  of  temperature  change  at  a  fixed  point  within  the  gas 
stream  is  insignificant  compared  to  the  change  of  temperature  of  a 
given  gas  element  as  it  moves  through  the  chamber,  it  is  possible 
to  rewrite  Eq.  49  as  follows 

*a*P  c  <T  )*4?““  =  -T4  U  <z.a,*,P'P  )  (51) 

8  P  g  02  8  1  88 


The  product  P 
to  yield 


z 

8 


may  be  replaced  by  the  term  for  mass  flow  W 


-T4  U.fz.a.-e.p'P  ) 
g  1  g  g 


(52) 


Due  to  the  fact  that  the  terms  c_  and  Uj  are  variable  in  T  and 
z  respectively,  this  equation  must  be  solved  numerically.  In  th% 
next  two  sections,  this  equation,  in  conjunction  with  other  equations 
and  information,  will  be  utilized  to  explore  the  effects  and  behavior 
of  burning  Composition  B  and  PBX  9404  under  pressures  of  a  few 
atmospheres.  This  will  be  followed  by  a  consideration  of  the  behavior 
of  burning  Composition  B  under  pressures  as  high  as  a  few  thousand 
atmospheres. 


(1)  Analysis  of  the  Effects  of  Burning  Composition  B  Under  Low 
Pressures 


This  section  will  be  concerned  with  the  determination  of  the 
constant  Pg  and  the  prediction  of  thermal  data  to  serve  as  a  check  on 
the  theory.  At  the  end  of  the  section,  an  analysis  will  be  conducted 
to  determine  the  minimum  size  of  vent  hole  which  will  prevent  an 
appreciable  pressure  buildup.  The  evaluation  of  pg  was  made  by 
adjusting  its  value  until  the  theoretical  relationship  between  heat  and 
mass  flow  was  satisfied.  To  accomplish  this,  it  was  necessary  to 
determine  the  temperature  distributions  of  the  gas  associated  with 
arbitrary  vaiuos  of  W,  Pg,  and  P_.  These  evaluations  were  made 
by  substituting  these  values  into  *q.  52.  In  each  case,  the  temperature 
of  the  gas  at  z  =  0.  was  determined  from  Eq.  46  and  Figures  2 
through  5.  The  radiant  flux  to  the  exploeive  was  calculated  by  sub¬ 
stituting  these  temperatures  into  Eq.  38  and  integrating  over  the 
several  gas  elements.  These  computations  were  performed  on  the 
UNIVAC  1105  computer  for  several  sets  of  parameters.  The  heat 
flux  which  is  conducted  to  the  explosive  at  the  assumed  pressure 
was  determined  from  the  data  shown  in  Figure  4  and  added  to  the 
radiant  flux  to  yield  the  total  heat  flux  to  the  explosive.  The  results 
were  then  checked  against  the  heat  flux  indicated  by  Eq.  45  for  the 
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assumed  mass  flow.  Appropriate  modifications  were  made  in  the 
parameters  until  the  relationship  between  heat  and  mass  flow  was 
closely  satisfied.  From  these  calculations!  several  sets  of  values 
of  the  parameters  W,  (3g  and  Pg  were  found  in  which  the  value  of 
(3g  was  2.83  x  10"4ft2  "cm*R/?t-lb  (pressure  in  units  of  lb/ft2  ).  This 
value  of  P'  appears  reasonable  since  it  indicates  that  the  gas  evolved 
by  the  burning  of  Composition  B  at  one  atmosphere  of  pressure  has 
an  opacity  which  is  comparable  to  that  of  flames  from  ordinary  hydro¬ 
carbon  fuels.  The  results  of  these  computations  are  shown  in 
Table  VII. 

The  fact  tnat  the  pressure  and  the  mass  flow  are  consistent 
with  test  results  will  be  discussed  at  the  end  of  this  section.  It 
may  be  observed  that  at  these  low  pressures  the  radiant  flux  is 
generally  several  times  larger  than  the  flux  due  to  conduction.  At 
pressures  of  the  order  of  tens  of  atmospheres,  the  heat  flux  due  to 
conduction  will  exceed  the  flux  due  to  radiation.  This  is  due  to  the 
fact  that  the  radiant  flux  gradually  levels  off  with  pressure  while  the 
conductive  flux  continues  to  increase. 

Before  applying  these  results,  it  is  desirable  to  check  various 
predictions  of  the  theory  against  experimental  results.  There  are 
three  types  of  thermal  data  which  may  be  used  as  a  heck  on  the 
theory.  These  are  the  temperature  distribution  ot  the  gas,  the 
irradiation  or  temperatures  associated  with  the  meters,  and  the 
temperature  distribution  along  the  length  of  the  chamber  at  various 
times. 

The  gas  temperatures  were  calculated  as  indicated  earlier  using 
Eq.  52  and  are  presented  in  Figure  29  for  the  several  burning 
velocities  noted  in  Table  VIII.  The  experimental  temperature  dis¬ 
tributions  corresponding  to  these  burning  velocities  were  found  by 
interpolating  between  the  data  shown  in  Figures  15  through  19. 

A  summary  of  the  theoretical  and  experimental  temperatures  is 
shown  in  Table  VIII,  It  may  be  observed  that  in  general  the  agree¬ 
ment  is  quite  good.  However,  it  also  should  be  noted  that  in  most 
cases  the  calculated  tempe ratures  exceed  the  experimentally  deter¬ 
mined  temperatures.  This  may  be  explained  by  the  difficulty  in 
experimentally  bringing  the  temperature  of  thermocouples  in  a 
partically  transparent,  low-velocity  gas  stream  up  to  the  temperature 
of  the  gas.  As  indicated  earlier  from  a  rough  calculation,  the 
experimental  values  could  be  low  by  as  much  as  50  to  200*F.  This 
would  account  for  the  pattern  of  the  discrepancies  shown  in  Table  VIII. 

The  first  step  in  the  evaluation  of  the  irradiation  to  the  meters 
was  the  substitution  of  the  gas  temperature  into  Eq.  38  and  replace¬ 
ment  of  the  term  £  by  the  distance  of  the  meter  (see  Figure  13) 
from  the  base  of  this  chamber.  The  irradiation  was  evaluated  by 
numerical  integration  of  the  expression  over  the  several  increments 
of  gas.  The  values  of  irradiation  were  then  substituted  into  Eq,  40, 
to  determine  the  temperatures  of  the  meters  for  the  several  meter 
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Table  VIII 

GAS  TEMPERATURES  {COMPOSITION  B) 


Burning  Velocity, 
ft/ hr 

Distance  from 
Explosive, 
ft 

Temperatures 

*F 

of  Gas, 

Experimental 

Calculated 

2.  5 

1.4 

1090 

1140 

1.6 

1040 

1050 

1.  8 

990 

980 

2.  0 

960 

910 

3.  5 

.  6 

1590 

1860 

.  8 

1450 

1650 

1.  0 

1320 

1500 

1.  2 

1210 

1350 

1.4 

1 1 10 

1250 

1.  6 

1040 

1160 

1.  8 

1000 

1080 

2.  0 

970 

1000 

5 

.  8 

1750 

1830 

1.  0 

1510 

1660 

i.  2 

1310 

1520 

1.4 

1170 

1410 

l.  6 

1080 

1300 

1.8 

1020 

1220 

2.  0 

990 

1140 

7.  5 

1.4 

1490 

1610 

1.6 

1350 

1500 

I.  8 

1250 

1400 

2.  0 

1180 

1320 

10 

1.  6 

1580 

1650 

1.8 

1450 

1550 

2.  0 

1370 

1460 

See  Figures  15  through  19. 
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locations  and  burning  velocities.  These  temperatures  are  presented 
in  Table  IX.  The  experimental  results  were  obtained  from  the  curves 
shown  in  Figure  20  to  23.  When  the  accuracy  of  the  meter  temper¬ 
atures  is  considered,  it  may  be  observed  that  the  agreement  is  very 
satisfactory. 


Another  means  of  checking  the  theory  is  to  compute  the  temperature 
distribution  along  the  length  of  the  chamber  (see  Figure  1)  and  com¬ 
pare  it  to  the  experimental  distribution  3hown  in  Figure  25.  This 
computation  is  much  more  complex  than  the  preceeding  ones  due  to 
the  fact  that  the  heat  flux  must  be  calculated  for  a  number  of  incre¬ 
ments-  of  the  chamber.  In  addition,  allowance  must  be  made  for  the 
fact  that  the  burning  front  is  moving,  causing  a  shift  of  the  gas 
temperatures  with  respect  to  the  chamber.  Allowance  must  also 
be  made  for  the  conduction  of  heat  along  the  length  of  steel  wall  of 
, the  chamber  during  and  after  the  burning. 

The  irradiance  of  the  elements  of  the  chamber  was  calculated 
by  first  modifying  Eq.  27.  This  equation  represents  the  rate  of 
heat  loss  from  a  particular  gas  element  to  the  total  curved  surface 
of  the  chamber.  For  purposes  of  the  computation,  it  is  necessary 
to  determine  the  total  heat  loss  from  all  the  gas  elements  to  a  parti¬ 
cular  ring-shaped  element  of  the  chamber  wall.  If  the  height  of  each 
chamber  element  and  the  thickness  of  each  gas  element  are  repre¬ 
sented  by  Az'  and  Az  (see  Figure  9),  then  by  manipulation  of  Eq.  27 

8P(z',a,P'P  ) 

Aq'^  =  T4  (z) - 8  Az Az '  (53) 

This  represents  the  radiation  which  is  emitted  by  a  gas  element 
and  is  incident  on  a  chamber  element  that  is  located  at  a  distance 
from  the  gas  element.  To  find  the  total  radiation  incident  on  the 
chamber  element  it  is  necessary  to  integrate  this  expression  over 
the  several  gas  elements.  If  b  and  £  are  used  to  represent  the 
distances  of  the  chamber  element  and  top  of  tire  chamber  from  the 
surf-ice  of  the  high-explosive  and  all  quantities  such  as  z',  z.  b 
and  £  are  taken  as  positive,  then 


(  b 
*>';  -  \  ■  f 


(54) 


/ 

J 


*  I  Vb, 


0  P  (z-b,  a,  (3 '  P  )  dz 
-  g  g 


Az' 
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Table  IX 

METER  TEMPERATURES  (COMPOSITION  B) 


Burning  Velocity, 
ft/hr 

Distance  from 
Explosive, 
ft 

Temperature  of  Meter, 

*F 

3.  5 

0.  42 

* 

Experimental 

1100 

Calculated 

990 

0.  82 

780 

730 

1.  22 

550 

580 

5.  0 

0.  42 

1200 

1160 

0.  82 

890 

890 

1.  22 

640 

720 

7.  5 

0.  42 

1480 

1390 

0.  82 

1040 

1090 

1,  22 

800 

900 

10.  0 

0.  82 

1220 

1240 

1.  22 

930 

1030 

See  Figures  20  through  23 
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This  represents  the  total  radiation  incident  on  a  chamber  ele¬ 
ment  that  has  a  thickness  Az1  and  is  located  a  distance  b  from  the 
high-explosive.  In  this  analysis  it  was  assumed  that  all  radiation  was 
absorbed  by  the  chamber  since  the  walls  were  blackened  by  the  soot. 
Also»  the  emission  of  radiation  from  the  surface  of  the  chamber  was 
neglected.  Convection  may  be  accounted  for  by  computing  the  heat 
transfer  coefficient,  h,  from  McAdams  9  and  applying  the  standard 
relationship  of 


Aq"  -  h(T  -  T  )  i  n aAz1  (55) 

g  w' 

Where  Tw  represents  the  wall  temperature. 

Due  to  the  approximate  nature  of  the  evaluation  of  the  heat 
transfer  coefficient,  the  value  of  h  was  considered  constant  along 
the  length  of  the  chamber  for  each  of  the  several  burning  rates.  The 
values  ranged  from  1  to  7  b/  ft2  hr*R. 

In  addition  to  the  radiant  and  convective  heat  fluxes  from  the  gas 
each  element  of  the  chamber  either  gains  or  loses  heat  due  to  heat 
conduction  along  the  length  of  the  steel  wall.  This  heat  flow  may  be 
represented  by  the  following  finite  difference  equation 


*q‘cond  =  -  Kstl  <  2a  +  d)*Z’  ( 

where 

Kgtl  -  thermal  conductivity  of  steel. 

As  a  result,  the  rate  of  change  in  temperature,  AT,  of  a  given 
element  of  the  chamber  may  be  determined  by  the  following  equation: 


PstlCstird(2a  +d,Az'^T  =^r  +*qc  +*qcond 


(57) 


where 


d 


Az* 

At 


specific  heat  of  steel,  B/lb*R 
density  of  steel,  lb/  ft3 
thickness  of  steel  wall,  ft 
radius  of  chamber,  ft 
height  of  element  of  chamber,  ft 
increment  of  time,  hr 


The  specific  heat  of  the  steel  was  varied  with  temperature  in 
accordance  with  the  values  shown  in  Table  VI  while  the  thermal  dif- 
fusivity,  i.  e. , 


vstl 


was  maintained  at  a  value  of  0.52ft/  hr2 


stl  stl 
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The  calculated  and  experimental  temperature  rises  of  the 
chamber  are  shown  in  Tables  X  and  XI  at  10  and  20  minutes  following 
the  start  of  burning  of  a  4-inch-long  cylinder  of  Composition  B  at 
the  base  of  the  chamber.  It  may  be  observed  that  the  agreement 
between  the  temperature  near  the  base  of  the  chamber  is  poor  while 
the  agreement  is  reasonably  good  at  distances  in  excess  of  about  0.  6 
feet  from  the  base.  The  low  experimental  temperatures  near  the  base 
of  the  chamber  are  attributed  to  the  presence  of  the  thermal  barrier 
effected  by  a  0.  050-inch-thick  coating  of  Adiprene  over  the  lower 
4  inches  of  the  chamber.  This  layer  of  Adiprene  was  bonded  both 
to  the  chamber  and  to  the  high- explosive  to  prevent  burning  along 
the  sides  of  the  Composition  B  cylinder. 


One  of  the  most  useful  predictions  evolved  by  this  theory  is 
the  minimum  size  of  the  vent  hole  which  will  prevent  a  significant 
pressure  buildup.  The  pressure  within  the  chamber  can  be  maintained 
at  relatively  low  levels  if  the  rate  of  mass  flow  through  the  vent  hole 
is  sufficient  to  counteract  the  flow  of  gas  evolved  by  the  burning  explo¬ 
sive.  To  determine  whether  or  not  this  condition  is  satisfied  requires 
an  evaluation  of  the  rates  at  which  gas  escapes  through  vent  holes  of 
various  sizes.  This  evaluation  may  be  accomplished  by  use  of  the 
following  two  equations10  describing  the  rate  of  mass  flow  W'  through 
a  vent  hole  which  presents  a  throat  area  of  At<  For  subsonic  flow 
through  the  throat 

r  '  •  *  2/k 


W' 


-  gJ 
R 

g 


for  sonic  flow 


W'  = 


A  P 
t  c 


AtPc«  J 


C 

T 


i  /,  k  +  1 

(*-)' -or 


1/2 


(58) 


k  +  1 


wha  re 


r  / 

v-rrr  ■) 

( 

k  +  1  , 

] 

1/2 


(59) 


W 

8 

J 

T 

c 

R 


rate  of  mass  flow  through  vent  hole,  lb/hr 

g 

gravitational  constant  =  4.  15  x  10°  ft/hr2 
mechanical  equivalent  of  heat  =  778  ft-lb/B 
temperature  of  gas  at  entrance  of  vent  hole,  *R 
gas  constant,  ft-lb/lb*R 
pressure  within  chamber,  lb/ft2 
pressure  at  throat,  lb/ft2 

ratio  of  "specific  heat"  at  constant  pressure  to  the 
"specific  heat"  at  constant  volume 

"specific  heat"  at  constant  pressure,  B/lb*R 
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Table  X 

TEMPERATURE  RISES  OF  CHAMBER  AT  10  MINUTES 
AFTER  START  OF  BURNING  COMPOSITION  B 


Distance  from 
Base  of  Chamber, 
ft 

Experimental 
Temperature  Rise, 

*F 

Calculated 
Temperature  Rise, 

•  f 

0 

405 

885 

0.  2 

495 

861 

0.  4 

590 

791 

0.  6 

625 

686 

0  8 

570 

567 

2.  0 

475 

455 

1,  2 

390 

359 

1.4 

305 

284 

1.6 

255 

232 

1.8 

215 

20  2 

2.  0 

195 

192 

RTD  TDR  63-3086 


Table  XI 

TEMPERATURE  RISES  OF  CHAMBER  AT  20  MINUTES 
AFTER  START  OF  BURNING  COMPOSITION  B 
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In  this  analysis,  the  throat  area  and  pressure  will  be  taken  as 
equal  to  the  cross-sectional  area  of  the  vent  hole  and  the  environmental 
pressure,  respectively.  It  may  be  observed  that  in  a  given  external 
environment,  the  flow  is  a  function  of  the  pressure  within  the  chamber 
and  the  temperature  of  the  gas  just  prior  to  entering  the  vent  hole. 

If  uniform  flow  through  the  chamber  is  assumed,  the  gas  temperature 
may  be  fixed  by  a  specification  of  the  pressure,  the  rate  of  burning, 
and  the  distance  between  the  high- explosive  and  the  vent  hole.  From 
this  information  the  gas  temperature  may  be  found  from  the  data 
shown  in  Table  VII  and  Figure  29.  In  order  to  determine  the  flow 
escaping  through  a  vent  hole,  it  is  first  necessary  to  determine 
whether  the  flow  is  subsonic  or  sonic.  This  may  be  determined  by 
the  following  equation*® 

k 


max 


(60) 


The  values  of  k  are  shown  as  a  function  of  gas  temperature  in 
Table  IV.  The  temperature  of  the  gas  entering  the  vent  hole  will 
range  from  1  370  to  5240*  R  for  Composition  B  with  vent-hole  distances 
between  0.  1  and  2.  0  feet  and  with  chamber  pressures  between  1  and 
5  atmospheres.  These  temperatures  will  yield  values  of  k  ranging 
from  1.  13  to  1.  36.  Substitution  of  these  values  into  Eq.  60  yields 
chamber  pressures  of  1.73  and  1.87  atmospheres.  Therefore,  for 
chamber  pressures  below  about  1.  73  atmospheres  the  flow  will  be 
subsonic  and  the  use  of  Eq.  58  is  indicated.  Chamber  pressures 
above  1. 87  atmospheres  will  produce  sonic  flow  at  the  throat  and 
the  use  of  Eq.  59  will  be  required. 

The  gas  temperatures  for  each  of  the  several  cases  were  ob¬ 
tained  from  the  data  shown  in  Table  VII  and  Figure  29.  The  depen¬ 
dence  of  R  and  k  on  temperature  may  be  found  in  Table  IV  while 
values  of  ’Tip"  may  be  determined  from  the  slopes  of  the  enthalpy 
temperature  curve  shown  in  Figure  5.  The  mass  flow  was  determined 
for  the  various  pressures,  distances,  and  vent  sizes  by  substitution 
of  appropriate  values  either  into  Eq.  58  if  the  flow  is  subsonic  or 
into  Eq.  59  is  the  flow  is  sonic.  The  results  are  shown  as  dashed 
curves  in  Figures  30  through  32.  The  solid  curve  in  each  figure 
refers  to  the  rate  at  which  gas  is  evolved  by  the  burning  Composition  B. 
By  comparing  the  figures,  it  may  be  observed  that  the  critical 
size  of  the  vent  hole  (for  which  the  inlet  flow  generated  by  the  burning 
explosive  is  always  equal  to  or  greater  than  the  exit  flow  passing 
through  the  vent)  increases  from  0.  28  to  0.  31  to  0.  35  inch  as  the  distance 
of  the  hole  from  the  explosive  decreases  from  2.  0  to  0.  5  to  0.  1  feet. 

With  the  larger  hole  sizes  shown  in  the  figures,  the  inlet  and  outlet 
flow  would  tend  to  stabilize  at  chamber  pressures  between  1.  1  to 
1.  2  atmospheres  and  result  in  mass  flows  of  about  30  to  40  ll/hr. 

This  range  of  mass  flows  corresponds  to  burning  velocites  between 
3.  3  and  4.  4  ft/hr  and  is  consistent  with  the  burning  velocities  of  about 
3.  5  ft/hr  that  were  encountered  in  tests  in  which  the  Composition  B 
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Fig.  30  INLET  AND  EXIT  MASS  FLOW  TO  CHAMBER 

(DISTANCE- COMPOSITION  B  AND  VENT  -  2FT) 
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burned  uneventfully.  In  addition,  the  chamber  pressures  are  in  agree¬ 
ment  with  the  values  indicated  above.  During  these  tests,  it  was 
frequently  observed  that  the  gas  and  smoke  escaped  in  a  pulsating 
manner.  This  may  be  attributed  to  pressure  oscillations  about  the 
pressure  indicated  by  the  first  intercept  of  the  dashed  curve  with  the 
solid  curve.  It  may  be  observed  that  if  the  pressure  is  lower  than 
this  value  the  flow  from  the  burning  Composition  B  will  exceed  the 
flow  escaping  through  the  vent  hole.  The  reverse  is  true  if  the 
pressure  exceeds  this  value.  In  either  case,  the  pressure  will  tend 
to  move  towards  the  value  indicated  above.  This  could  cause  oscil¬ 
lations  of  the  pressure  and  burning  rate  which  would  explain  the  puffs 
of  escaping  gas  and  smoke.  The  amplitude  of  the  pressure  oscillations 
must  have  been  less  than  about  0.  1  atmospheres  since  no  oscillations 
were  detected  by  a  catenary  type  pressure  gage  whose  sensitivity  was 
about  0.  1  atmosphere.  The  experimental  data  shown  in  Figure  24 
also  provide  a  check  on  the  results.  Those  data  indicate  that  the 
critical  diameter  of  the  vent  hole  associated  with  the  2  foot  chamber 
shown  in  Figure  13  lies  between  0.  222  and  0.  304  inch,  and  is  in  good 
agreement  with  the  value  of  0.  28  inch  shown  in  Figure  30. 

As  noted  previously  the  pressure  may  be  maintained  at  tolerably 
low  levels  by  use  of  vent  holes  whose  diameters  exceed  the  critical 
values  of  0.  28  to  0.  35  inch  associated  with  vent-explosive  distances 
of  0.  1  to  2.  0  feet,  respectively.  The  cross-sectional  areas  of  these 
holes  are  equivalent  to  0.  49  to  0.  77  percent  of  that  of  the  burning 
surface.  It  may  be  observed  that  the  critical  area  of  the  vent  hole 
will  increase  according  to  the  proximity  of  the  vent  hole  to  the  explo¬ 
sive.  Except  for  minor  differences,  these  data  may  be  applied  to 
systems  employing  Cyclotol  type  explosives  for  which  the  burning 
rate  and  gas  temperatures  are  nearly  the  same  as  those  for  Com¬ 
position  B.  If  these  data  are  to  be  applied  to  larger  systems,  such 
as  actual  weapon  systems,  care  should  be  taken  to  provide  adequate 
coverage  of  the  pressure  hull  with  the  vent  holes  to  ensure  that  no 
significant  pressure  gradients  develop  as  a  result  of  internal  resis¬ 
tances  to  gas  flow. 

(2)  Preliminary  Analysis  of  the  Effects  of  Burning  PBX  9404  Under 

Low  Pressures 

As  was  noted  in  Section  3,  PBX  9404  burns  approximately 
twice  as  fast  as  Composition  B  at  atmospheric  pressures.  The 
faster  burning  rate  is  attributed^  to  the  fact  that  PBX  9404  reacts  more 
completely  than  Composition  B.  This  results  in  higher  gas  temperatures 
and  hence,  higher  radiant  fluxes  to  the  explosive. 

Unfortunately,  due  to  the  lack  of  strand  burn  data  for  either 
PBX  9404  or  HMX  which  is  its  major  constituent,  it  was  not  possible 
to  accurately  determine  the  heat  flux  to  this  explosive  due  to  con¬ 
duction.  For  this  reason,  the  study  was  confined  to  a  preliminary 
examination  of  the  heat  flux  based  on  the  temperatures  of  the  gas  and  of  the 
chamber  shown  in  Figure  1.  This  analysis  indicates  that  the  value 
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of  (3A  is  close  to  0.47  x  10”  ft2  cm*R/ft-lbi  while  the  rate  of 
conductive  heating  at  a  pressure  of  1  atmosphere  is  in  the  neighbor¬ 
hood  of  75#  000  B/hr-ft2  .  The  fact  that  the  value  of  is  appreciably 
less  than  the  value  for  Composition  B,  i.  e. ,  2.83  x  1(T4.  indicates 
that  the  flames  from  PBX  9404  are  much  cleaner  and  transparent 
than  those  of  Composition  B.  This  is  consistent  with  the  observation 
that  the  burning  of  Composition  B  results  in  considerably  more  soot 


than  the  burning  of  PBX  9404, 
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4.  BURNING  OF  COMPOSITION  B  AND  PBX  9404  UNDER  HIGH 

PRESSURES 

The  experimental  studies  were  performed  by  burning  Com¬ 
position  B  and  PBX  9404  under  pressures  ranging  from  one  to 
several  thousand  atmospheres.  In  addition,  an  analysis  was  conducted 
to  predict  the  pressure  using  the  theoretical  and  experimental  results 
discussed  previously.  The  purpose  of  this  work  is  two-fold:  first 
to  see  if  it  is  possible  to  predict  the  transient  pressure  rises  during 
normal  burning;  and  second,  to  search  for  analytical  criteria  which 
could  be  used  to  predict  the  type  or  order  of  explosion  which  may  occur. 
The  first  objective  was  successfully  accomplished  and  the  second 
objective  was  only  partially  successful. 

The  remainder  of  this  section  will  discuss  the  results  of  the 
experimental  work  and  will  present  an  analytical  means  of  pre¬ 
dicting  the  pressure  history  during  normal  burning. 

a.  Experimental  Results  from  High  Pressure  Burning  of 

Composition  B 

This  section  will  cover  the  burning  of  Composition  B  under 
pressures  up  to  several  hundred  atmospheres.  The  pressures  were 
obtained  by  burning  the  high-explosive  in  tubular  chambers  which 
were  sealed  to  prevent  gas  leakage  and  which  were  strengthened  to 
maintain  pressures  of  the  order  of  100  atmospheres.  A  sketch  of 
the  chamber  assembly  is  shown  in  Figure  33.  Photographs  of  one 
of  the  steel  chambers  are  shown  in  Figures  34  and  35.  The  space 
of  about  0.  05  inches  between  the  Composition  B  and  the  wall  of  the 
chamber  was  sealed  with  Adiprene  to  prevent  ignition  of  the  sides 
of  the  Compos itioi.  B  cylinder  by  exposure  to  the  hot  gas.  The  dia¬ 
meter  of  the  rods  and  the  thickness  of  the  cover  plates  were  increased 
above  the  values  used  in  the  low-pressure  chamber  shown  in  Figure 
13  to  cause  very  violent  explosions  or  low  order  detonations.  During 
assembly  of  the  chamber,  the  rods  were  elongated  to  within  70  per¬ 
cent  of  their  yield  strength  in  order  to  retain  the  cover  plates  for  as 
long  is  practical.  This  was  accomplished  by  measuring  the  elongation 
of  the  rods  with  a  dial  gage  as  torque  was  applied  to  the  nuts.  The 
upper  face  of  the  explosive  was  heated  and  ignited  using  the  same 
technique  illustrated  previously  for  the  low-pressure  tests.  This 
was  assomplished  by  passing  a  current  through  the  122-watt  heater 
shown  in  Figure  14  for  a  period  of  10  minutes.  This  initiated  melting 
of  the  surface  of  the  Composition  B.  Ignition  was  accomplished  by 
passing  a  current  of  10  amperes  through  a  0.56-ohm  wire  which  was 
located  in  a  shallow  groove  across  the  face  of  the  Composition  B.  An 
average  delay  of  5.  4  seconds  was  noted  between  the  delivery  of  this 
current  and  the  start  of  burning  (detected  by  an  abrupt  decrease  of  the 
current  through  the  igniter  wire).  Various  types  of  ionization  probes 
were  used  to  measure  the  passage  of  the  flame  front  across  various 
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Fig.  33  CROSS  SECTION  OF  "CLOSED"  TEST  CHAMBER 
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Fig.  34  SIDE  VIEW  OF  "CLOSED"  CHAMBER  FOR 
HIGH-PRESSURE  BURNING 
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Fig.  35  BOTTOM  VIEW  OF  "CLOSED"  CHAMBER  FOR 
HIGH-PRESSURE  BURNING 


I 

j 
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points  along  the  axis  of  the  explosive.  Sketches  of  two  of  the  most 
commonly  used  probes  are  shown  in  Figures  36  and  37.  The  probes 
were  inserted  into  the  Composition  B  through  small  holes  in  the 
chamber  wall.  The  technique  of  positioning  the  probes  and  sealing 
the  holes  is  illustrated  in  Figure  38.  This  figure  also  contains  a 
sketch  which  shows  the  method  used  to  mount  a  Kistler  pressure 
gage  in  the  wall  of  the  chamber. 

During  the  course  of  this  work,  a  total  of  30  tests  were  con¬ 
ducted  in  which  the  distance  between  the  explosive  and  the  top  of 
the  chamber  was  varied  from  l/ 4  f.o  31  inches  while  the  minor  dia¬ 
meter  of  the  rods  was  varied  from  0.  1  9  to  0. 85  inches.  Two  types 
of  burning  were  noted  which  were  distinguished  from  each  other  by  a 
marked  difference  in  the  rates  of  pressure  rise  within  the  chamber. 
Illustrations  of  this  may  be  seen  in  Figures  39  through  41.  These 
figures  show  the  pressure  histories  of  three  tests  in  which  the  dis¬ 
tance  between  the  high-explosive  and  the  lower  surface  of  the  top- 
cover  plate  was  varied.  The  portions  of  the  curves  between  14.  7 
and  about  65  psia  were  omitted  due  to  the  appreciable  times  required 
for  this  transition  as  compared  to  the  time  scales  shown  in  the 
figures.  In  each  case,  it  may  be  observed  that  the  pressure  built 
up  relatively  slowly  to  about  1000  psia  and  then  rose  very  abruptly 
as  evidenced  by  the  discontinuity  of  the  slope  of  the  pressure-time 
curve.  A  similar  ano  naly  was  noted  by  WachtsU^  with  smaller 
1  inch  diameter;  Composition  R  samples  at  pressure-  of  4000- 
5000  psia.  The  portion  of  the  pressure-time  curve  up  to  about 
1000  psia  is  consistent  with  what  would  be  expected  if  the  burning  is 
confined  to  the  face  of  the  Compos iticr.  B  cylinder.  This  will  he 
illustrated  in  the  next  sections  The  very  abrupt  rise  in  pressure 
is  attributed  to  an  increase  of  the  burning  area  caused  by  cracking 
of  the  explosive. 

There  are  three  pieces  of  evidence  which  indicate  that  the 
Composition  B  cracked.  The  most  conclusive  evidence  is  illustrated 
by  a  photograph  of  a  number  of  pieces  of  Composition  B  shown  in 
Figure  42.  These  pieces  of  explosive  were  recovered  from  a  test 
in  w.-’ch  the  top  cover  plate  was  ejected  and  caused  the  burning  to  be 
extinguished  due  to  the  rapid  decrease  of  gas  pressure  and  temperature. 
The  porous  surfaces  arc  those  surfaces  which  were  formed  during 
burning  and  in  the  process  of  very  vigorous  reaction.  It  may  be 
observed  that  this  porous  condition  exists  on  irregular  surfaces  of 
practically  all  the  pieces  o'  Composition  E  and  could  happen  only  if 
the  high-explosive  cracked  and  exposed  these  surfaces  to  hot,  high- 
pressure  gas  prior  to  ejection  of  the  cover  plate.  The  second  piece 
of  evidence  that  indicates  the  Composition  B  cracked  is  the  random 
sequence  in  which  the  probes  fired.  This  happened  with  each  of  the 
several  types  of  probes.  Extensive  cracking  of  the  high- explosive 
would  result  in  random  firing  of  the  probes  and  appears  to  be  the 
most  likely  cause  of  the  irregular  firing  sequence.  The  last  evidence 
which  supports  the  theory  of  cracking  is  the  apparent  susceptibility 
of  Composition  B  to  thermal  shock.  The  Composition  B  was  stored 


79 


RTD  TDR  63-3086 


Fig.  37  TYPE-B  ION  PROBE  AND  ASSOCIATED  CIRCUIT 


81 


RTD  TDR  63-3086 


PICKUP  (601)  MOUNTED  INTO  CHAMBER  WALL 


ASSEMBLY  OF  PRESSURE  SEAL  FOR 


THERMOCOUPLE  LEADS 


Fig.  38  METHOD  OF  INSERTING  THERMOCOUPLES  AND 
PRESSURE  PROBE  THROUGH  CHAMBER  WALL 
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PRESSURE  TIME  HISTORY  (COMPOSITION  B,  i  =  1.  75  inches) 
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in  an  igloo  which  often  reached  temperatures  of  0*F  or  less.  Upon 
exposure  of  the  high- explosive  to  normal  room  temperatures,  it  com¬ 
menced  to  emit  audible  cracking  noises  although  no  visual  changes 
were  discernable.  When  sawed,  the  Composition  B  cracked  into  a 
number  of  pieces.  This  gives  an  indication  of  the  susceptibility  of  the 
explosive  to  thermal  shock.  During  the  burning,  the  solid  explosive 
beneath  the  layer  of  foam  is  subjected  to  very  appreciable  temperature 
gradients.  An  appreciation  of  the  magnitude  of  these  temperature 
gradients  may  be  gained  by  a  consideration  of  their  magnitude  during 
constant- velocity  burning.  Under  this  condition,  the  temperature,  T, 
of  the  solid  Composition  B  may  be  described  by  ^ 


T 


where 


u 

a 

x 


(T  -  T  )  exp  -(-£-  x)  +  T 
m  o  n  o 

melt  temperature  of  high-explosive,  °F 
initial  temperature  of  high-explosive,  *F 
burning  velocity,  ft/hr 

thermal  diffusivity  of  solid  high-explosive,  ft2 /hr 
distance  from  melt  surface,  ft. 


(61) 


The  temperature  gradient  at  the  surface  may  be  found  by  differ¬ 
entiating  this  expression  with  respect  to  x  and  setting  x  -  0  to  yield 


dT 

dx 


(62) 


From  Table  A-l,  Tm  and  a  are  approximately  equal  to  176*F 
and  .  0044  ft2 /hr  where  the  value  of  the  thermal  conductivity  of  Cyclotol 
70/30  was  used  in  the  evaluation  of  a.  Therefore,  if  the  initial  temper¬ 
ature  of  the  explosive  is  70* F,  the  temperature  gradient  at  the  surface 
is 


dT  _ 
*dx"  " 


-2.4  x  10 


(63) 


Even  for  the  lowest  burning  velocities  of  2.  5  ft/hr,  the  gradient 
is  substantial,  i,  e. ,  -6  x  104  *F/ft.  The  magnitude  of  the  temperature 

gradient  will  increase  in  a  linear  fashion  with  the  velocity  of  burning 
and  probably  contributes  to  the  cracking.  An  exhaustive  consideration 
of  the  problem  of  cracking  is  outside  the  scope  of  this  investigation. 
Future  studies  are  needed  in  this  area. 


The  effect  of  the  pressure  rise  on  the  chamber  varied  according 
to  the  strength  of  the  weakest  portion  of  the  chamber,  i.  e. ,  the  8 
ordinary  steel  rods  holding  the  cover  plates.  These  tests  were  con¬ 
ducted  utilizing  rods  with  the  following  threads  and  minor  diameters: 
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Thread  Size 
l/4  x  20 
5/16  x  18 
8/8  x  16 
7/16  x  14 
1/2  x  13 
5/8  x  11 
1  x  8 


Minor  Diameter 
.  189  inches 
,  244  inches 
■  298  inches 
.  3  5  C  inches 
.  406  inches 
.  514  inches 
.  84 7  inches 


With  small  threaded  rods  having  a  minor  diameter  of  0.  189  inches, 
the  pressure  reached  levels  of  several  hundred  psia  and  raised  the  top 
cover  plate  by  either  elongating  or  breaking  the  rods.  The  smallest 
rod  which  caused  destruction  of  the  chamber  had  a  minor  diameter  of 
0.  244  inch.  This  would  indicate  that  the  critical  value  of  the  minor 
diameter  (which  may  or  may  not  cause  destruction  of  the  chamber)  lies 
somewhere  between  0.  189  and  0.  244  inch.  The  static  pressure  which 
will  cause  separation  of  the  cover  plate  from  the  steel  tube  for  these 
two  rod  sizes  is  only  slightly  less  than  the  dynamic  pressure  at  separation 
and  may  be  determined  by  the  following  equation 


P  =  static  pressure  to  cause  separation  of  cover  plate  from 
steel  tube,  psia 

Y  =  preload  fraction  of  yield  stress  =  .  7  (dimensionless) 

<r  =  yield  stress  of  bolts  =  75,  000  psia 

y 

N  “  number  of  bolts  =  8 

=  minor  radius  of  bolts,  in. 

A  =  area  of  tube  face  =  7.  07  in.  2 

c 

A  =  area  of  cover  pluie  subjected  to  gas  pressure  =  12.57  in.  2 

Substitution  of  these  quantities  and  the  radius  of  the  critical  sizes 
of  bolts  into  this  equation  indicates  that  the  critical  pressure  necessary 
to  initiate  the  destruction  of  the  chamber  lies  somewhere  between  960 
and  1650  psia.  This  conclusion  is  reasonably  consistent  with  the  evidence 
which  indicates  that  normal  Composition  B  cracks  during  burning  upon 
exposure  to  pressures  in  the  vicinity  of  1000  psia. 

With  the  large  threaded  rods  that  had  a  minor  diameter  of  0.51  inch 
or  greater,  the  chamber  was  always  destroyed  by  the  reaction.  The 
damage  varied  from  test  to  test  with  rods  having  intermediate  dia¬ 
meters.  In  approximately  l/4  of  the  tests  utilizing  rods  with  minor 
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diameters  of  0,  24  to  0.  51  inch,  the  chamber  was  destroyed.  In  the 
remaining  testa,  the  damage  was  principally  confined  to  the  rods, 
even  though  the  pressure  exceeded  1000  psia  or  68  atmospheres  in 
each  of  these  tests.  This  anomaly  may  be  explained  by  variations 
in  the  extent  of  cracking  of  the  explosive  which  resulted  in  significant 
variations  in  the  increase  of  burning  area  and/or  the  development  of 
high  pressures  or  shocks  between  crevices.  This  of  course,  would 
cause  variations  in  the  violence  of  the  explosions.  Photographs  of  the 
remains  of  two  chambers  are  shown  in  Figures  43  and  44.  Figure 
43  represents  a  plastic  failure  of  the  chamber  in  which  the  walls 
were  thinned  out  and  ruptured  along  lines  extending  along  the  length 
of  the  chamber.  The  damage  was  typical  of  the  lower  order  explosions. 
The  effects  of  a  high  order  explosion  are  shown  in  Figure  44.  This 
debris  represents  the  remains  of  the  chamber  shown  in  Figures  34 
and  35.  The  inner  surface  c£  a  rumber  of  the  pieces  of  the  chamber 
were  wrinkled  as  shown  by  the  photographs  at  the  bottom  of  Figure  45. 
This  gives  an  indication  of  the  violence  of  the  explosion.  The  pres¬ 
sure  within  the  chambers  that  were  destroyed  must  have  been  con¬ 
siderably  higher  than  the  1300  atmospheres  required  to  cause  plastic 
failure  of  the  ordinary  carbon  steel  tube  under  static  loads. 

The  relatively  slow  response  time,  i.  e.»  2.8  milliseconds,  of  the 
galvonometers  in  the  CEC  recorder  did  not  resolve  the  high  pressure 
por!  .on  of  the  curves,  which  included  these  peak  pressures. 

b.  Experimental  Results  from  High  Pressure  Burning  of  PBX  9404 

Three  high  pressure  teste  were  conducted  using  PBX  9404. 

Except  for  the  omission  of  the  heater  which  was  not  deemed  necessary 
to  initiate  a  vigorous  burning,  the  tests  were  identical  to  those  with 
Composition  B.  The  chamher  was  12  inches  long  in  each  of  the  three 
tests,  and  the  rods  had  1/4-20  and  3/8-16  threads.  Within  times 
of  the  order  cf  milliseconds  following  the  ignition  of  5  inch  long 
samples  of  PBX  9404,  the  pressure  built  up  and  destroyed  the  chamber 
and  stand,  A  photograph  of  the  remains  of  one  of  the  chambers  is 
shown  in  Figure  46.  It  may  be  observed  that  all  components  of  the 
chamber  were  fractured  and  reduced  to  many  more  pieces  than  in 
the  Composition  B  tests  (see  Figures  43  and  44).  Plastic  flow  along 
the  inner  surface  of  the  chamber  was  noted  on  a  number  of  the  pieces. 
These  observations  indicate  that  the  reaction  was  a  detonation. 

Due  to  the  wide  variation  in  the  rise  times  of  the  pressures, 
i.  e.  ,  by  a  factor  of  at  least  10,  it  is  suspected  that  these  samples 
cracked.  This  would  account  for  large  differences  in  burning  since 
it  would  not  be  expected  that  each  sample  would  crack  in  the  same 
fashion.  Unfortunately,  due  to  the  very  rapid  rise  times  relative  to 
the  response  time  of  the  CEC  recorder  it  was  not  possible  to  discern 
if  and  when  the  slope  of  the  pressure  time  curve  exhibited  a  discontinuity 
which  would  indicate  that  the  PBX  underwent  extensive  cracking.  It 
is  recommended  that  any  future  studies  involving  PBX  type  explosives 
utilize  much  faster  response  instruments  such  as  oscilloscopes  to 
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Fig.  43  PLASTIC  FAILURE  OF  CHAMBER  (COMPOSITIONS) 
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EFFECT  OF  EXPLOSION  ON  INNER  SURFACES 
OF  CHAMBER  (COMPOSITION  B) 
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detect  any  occurence  of  cracking.  Elaborate  triggering  circuits  will 
be  required  to  ensure  that  the  pressure  is  recorded  to  allow  for  the  large 
differences  in  the  time  scale. 

c.  Theoretical  Analysis  of  the  Effects  of  Transient  Burning  of  High- 

Explosives  (High  Pressure) 

This  section  describes  the  development  of  an  analytic  description 
of  the  pressure  rise  within  the  closed  chamber  shown  in  Figure  34. 

The  work  was  programmed  for  the  UNIVAC  1105  computer  in  a  man¬ 
ner  that  allows  application  to  a  wide  variety  of  problems  of  this  type. 

The  temperature  distribution,  Tfc(x),  along  the  axis  of  the 
explosive  at  the  time  current  is  passed  through  the  igniter  wire  will 
be  represented  by 

Tf„M  *  Tm»*P  +  To  <65> 

where 

Tm  is  the  temperature  of  the  high-explosive  at  x  =  0.  *R 

Tq  is  the  initial  temperature  of  the  high- explosive  prior  to  any 
heating,  *R 

Cj  is  an  arbitrary  constant,  l/ft 

At  the  time  current  was  passed  through  the  igniter  wire,  the 
surface  was  in  the  process  of  melting.  An  average  time  of  5.4 
seconds  was  required  to  initiate  a  burning  reaction.  The  heat  from 
the  igniter  would  be  dissipated  in  depth  as  well  as  along  a  portion  of  the 
molten  film  on  the  face  of  the  Composition  B,  In  order  to  maintain 
a  simple  geometric  description  of  the  problem,  the  heat  will  be  as¬ 
sumed  as  dissipated  uniformly  in  a  thin  layer  across  the  face  of  the 
Composition  B,  The  thickness  of  the  layer  will  be  taken  equal  to  the 
thickness  of  the  foam  layer  during  the  burning  reaction  while  the 
width  of  the  element  will  he  adjusted  to  yield  ignition  times  consistent 
with  the  average  delay  of  5.  4  seconds  noted  above.  The  area  of  the 
upper  surface  of  the  element  will  be  represented  by  Aj  while  the 
remaining  area  of  the  face  of  the  Composition  B  will  betaken  as  A^. 

The  total  area  of  the  face,  i.  e. ,  Aj  +  A2  will  be  represented  by  A. 

The  gas  in  the  chamber  at  the  time  current  is  passed  through  the 
igniter  wire  will  be  considered  aB  having  the  same  properties  as  the 
gas  evolved  by  the  explosive.  The  total  mass  M  of  the  ga3  in  the 
chamber  may  be  represented  as  follows: 

M  =  (AjXj  +  A2x£  +  Axq)Pe  (66) 

where 

Xj  =  distance  traversed  by  the  burning  front  through  the 
element  of  the  high-explosive  containing  the  igniter 
wire,  ft 

=  distance  traversed  by  the  burning  front  through  the 
remainder  of  the  high-explosive,  ft 
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x.  =  depth  of  high- explosive  presenting  the  same  mass  as 

gas  originally  in  the  chamber,  ft 
Pj,  =  density  of  solid  high- explosive,  lb/ft1 

Since  the  volume  of  the  gas  is  essentially  constant,  the  rate  of 
energy  change  of  the  gas  may  be  represented  by 


.  /  H(T  )  \  cbc,  1 

Ve  -ar-  +  HiITn»j  (67) 

+  qelec  +  A2PE  ~W  /^2*X2J  +  H2(TflJJ  ■AlqcrA2qc2‘l,ri 

whe  re 

H(T  )  represents  the  total  equilibrium  enthalpy  of  the  gas,  B/lb 

k  =  ratio  of  the  "specific  heat"  at  constant  pressure  to  that 
at  constant  volume,  dimensionless 

q.(x.)  =  heat  released  by  the  reaction  of  material  evolved  from 
■  the  ith  element  of  explosive  (i  =  1,  2),  B/lb 


Hi(Tfi)= 


^elec 

qci 

qif  = 


enthalpy  of  explosive  of  the  ith  element  (  i  =  1,-2)  just 
prior  to  reaction,  B/lb 

rate  at  which  gas  is  heated  by  the  heater,  B/ft2 

rate  of  heat  extraction  from  gas  due  to  conduction  to 
the  ith  element  of  the  explosive  (i  =  1,  2),  B/ft2  hr 

rate  of  heat  extraction  from  gas  due  to  radiant  losses 
to  the  chamber  and  explosive,  B/ft2  hr 


If  the  gas  is  ideal,  the  equation  of  state  may  be  expressed  as 


P  (  £+x)  r  (X  +  x Jp_R  T  (68) 

8  °  E  g  8 

whe  re 

£  =  length  of  gas  volume,  ft 

From  Eq,  48  the  rate  of  radiant  heat  loss,  q  £»  from  the  gas  is 
represented  by 

£ 

q  „  =  |  T*  (z)  U  <z,a,  £,P'P  )dz  (69) 

rx  g  l  8  8 

o 

while  the  rate  of  heat  flow,  q  .  to  the  ith  element  of  the  explosive  is  given 
by  gl 
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q8i 


£. 

J  T*g(«)S1(«.a,l.P^Pg)dz  +  A.qci  +  . 


(70) 


where 

^  eleci  1 

^  elec.  2 


rate  of  heat  deposition  into  element  1  from  the 
heater  and  from  the  igniter,  B/lb 

rate  of  heat  deposition  into  element  2  from  the 
heater,  B/lb 


This  heat  flow  will  be  equal  to  the  rate  at  which  heat  is  accumulated 
in  the  ith  element  of  explosive  plus  the  different  between  the  rates  at 
which  heat  is  carried  into  and  out  of  the  layer  by  the  material  entering 
and  leaving  the  element.  This  may  be  expressed  as  follows: 


*gx 


whc  rc 


Ve 


dx. 

i 

nr 


H(Tfi)  -  H(Tfo.) 


+  M. 


dH(Tfi) 

"^eT 


dTf. 

fi 

ST 


temperature  of  foam  in  ith  element  during  burning,  *R 

temperature  of  solid  HE  just  prior  to  entering  foam  layer 
of  ith  element 


This  equation  assumes  that  the  amount  of  high-explosive  material 
in  the  foam  layer  remains  constant  and  that  the  solid  HE  under  the 
foam  layer  remains  at  its  initial  temperature  until  it  enters  the  foam 
layer.  Any  heat  abstracted  from  the  foam  by  evaporization  is  neglected. 


From  Eq.  4  the  rate  at  which  mass  is  evolved  by  the  ith  element 
of  the  explosive  (i  =  1,  2)  is 

M.  =  M.Z  exp  -(E/RgQTf.)  (71) 

Thia  completes  the  specification  of  the  several  phenomena.  All 
the  quantities  are  well  defined  except  for  the  heat  flow  to  the  explosive 
by  conduction,  i.  e. ,  q0j.  During  the  time  the  explosive  is  burning, 
the  values  of  qcj  are  defined  by  the  data  shown  in  Figure  4.  The 
problem  is  to  specify  when  these  data  are  applicable,  i.  e. ,  when 
burning  commences.  The  temperature  of  the  foam  layer  during  steady 
burning  at  a  pressure  of  1  atmosphere  is  approximately  540*F.  If 
one  uses  this  temperature  as  a  criterion  of  ignition  it  will  be  found  that 
the  pressure  at  which  the  pressure  commences  to  rise  rapidly  is 
much  higher  than  the  value  found  experimentally,  i,  e. ,  by  a  factor 
of  about  2.  This  would  indicate  that  burning  started  at  a  lower  foam 
temperature.  By  trial  and  error  methods  it  was  found  that  a  good 
criterion  to  prescribe  the  start  of  burning  is  the  time  at  which  the 
foam  temperature  reaches  510*F.  At  this  time  the  theoretical  pres¬ 
sure  is  consistent  with  the  experimental  values  at  which  the  pressure 
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commences  to  rise  rapidly.  Furthermore,  this  foam  temperature  is 
consistent  with  the  surface  temperatures  of  Composition  B  and  Cyclotol 
type  explosives  that  are  generally  encountered  In  model  and  full- 
scale  tests  just  prior  to  a  rapid  temperature  rise. 


These  results  were  programmed  for  the  UNIVAC  1105  computer  to 
determine  the  pressure  and  temperature  histories  of  the  gast  the 
temperature  of  the  layers  of  explosives(  and  the  burning  velocity  of 
the  Composition  B  in  the  chamber  shown  in  Figure  34.  In  addition, 
the  temperature  histories  of  the  chamber  wall  and  the  top  cover 
plate  were  computed.  The  parameters  used  in  these  computations 
are  as  follows 


M 

M 

R 


1 

2 


8 


47,  500  cal/ mole 
universal  gas  constant 

frequency  factor 


=  1.  987  cal/mole*K 

=  1.  104  cal/mole*R 

=  1018,5/sec 

=  1.  14  x  102  2  /hr 


0.  0238  lb 
0.  0491  lb 
73  ft- lb/lb  *R 

0.6  atmos.  lb2  *R  ft/ {ft- lb )2  cm  or  2.83  x  10*4lb2  *R/ 
ft2  ft- lb  cm 


k 

^elec 

qelec 

_  II 

^  elec 


1.31  (dimensionless) 

heat  absorbed  by  gas  from  heater  =  83B/hr 

rate  of  heating  of  high-explosive  from  igniter  wire 
=  209  B/hr 

rate  of  heating  of  high- explosive  by  heater  =  1.9x10 
ft2  hr 


B/ 


Three  calculations  were  conducted  utilizing  the  chamber  dimensions 
asso.  fated  with  the  three  tests  indicated  by  the  data  shown  in  Figures  39 
through  41.  In  these  tests,  the  distances  between  the  Composition  B 
and  the  lower  surface  of  the  top  cover  plate  (see  Figure  33)  were 
1.  75,  8.  0  and  17.  5  inches.  The  calculated  and  experimental  pres¬ 
sures  are  shown  in  Tables  XII  through  XIV.  It  may  be  observed  that 
the  agreement  between  the  experimental  and  calculated  pressures 
is  reasonably  good  up  to  about  1000  psia.  Beyond  this  point  the  experi¬ 
mental  pressures  are  considerably  higher  than  the  calculated  values. 
This  is  attributed  to  the  increase  of  the  burning  area  caused  by  the 
cracking  of  the  explosive.  While  these  data  are  consistent,  it  should 
not  be  expected  that  all  Composition  B  samples  will  crack  at  about 
1000  psia.  The  fact  tiiat  the  1-1  l/4  inch  diameter  samples  of  Com¬ 
position  B  cracked  in  the  Picatinny  Arsenal  tests^  at  4000  -  5000  psia 
suggests  that  size,  as  well  as  other  parameters,  such  as  severe 
cycling  of  the  environmental  temperature  and  rate  of  pressure  rise, 
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may  play  a  significant  role  in  affecting  the  cracking.  For  these  reasons, 
the  criterion  of  cracking,  i.  e.»  1000  psia,  is  only  preliminary  and 
should  be  treated  accordingly  until  the  cause  is  thoroughly  explored. 

Tables  XII  through  XIV  also  show  the  calculated  temperatures  of 
the  foam,  gas  and  chamber.  It  may  be  observed  that  at  the  lower 
pressures  there  is  a  substantial  difference  between  the  temperature 
of  the  foam  in  -lose  proximity  to  the  igniter  wire  and  the  temperature 
of  the  remainder  of  the  foam.  This  would  indicate  that  the  burning 
at  the  lower  pressures  was  essentially  limited  to  the  high-explo&ive 
surrounding  the  igniter  wire.  At  later  times,  the  computation  indicates 
that  the  complete  face  of  the  Composition  B  burned  at  the  same  rate. 

The  amount  of  material  that  is  burnt  away  during  the  transition  was 
quite  small,  i.  e..  of  the  order  to  0.  1  inch.  However,  since  much  of 
this  material  was  evolved  in  close  proximity  to  the  igniter  wire,  there 
is  a  possibility  that  a  groove  may  have  been  formed  across  the  face 
of  the  Composition  B  which  could  initiate  the  cracking.  This  would 
explain  the  difference  between  the  pressures  of  1000  psia  at  which 
the  Composition  B  cracked  in  these  studies  and  the  pressure  range 
of  4000  -  5000  psia  found  at  Picatinny  Arsenal.  This  point,  however, 
remains  to  be  resolved. 

The  calculations  also  revealed  that  the  average  temperature  rise 
of  the  chamber  is  relatively  small  in  spite  of  the  fact  that  the  heat 
fluxes  are  very  high.  This  is  due  to  the  small  times  for  heat  trans¬ 
fer.  The  fact  that  the  temperature  rises  of  the  chamber  are  only 
4*  to  16*F  at  the  time  of  cracking  for  the  different  tests  leads  one 
to  suspect  that  the  thermal  expansion  of  the  chamber  was  not  a  primary 
cause  of  the  high-explosive  cracking.  It  remains  to  be  seer,  whether 
or  not  the  expansion  of  the  chamber  by  the  pressure  played  any  role 
in  the  cracking  of  the  Composition  B. 
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5.  SUMMARY  AND  CONCLUSIONS 

This  report  indicates  it  is  feasible  to  make  reasonably  accurate 
predictions  of  the  buildup  of  gas  pressure  within  systems  containing 
burning  high-explosives  at  least  until  the  hig'h- explosive  cracks.  For 
this  purpose,  computer  programs  have  been  developed  to  determine 
the  pressure  within  closed  and  vented  systems.  For  these  computations, 
it  is  necessary  to  know  the  enthalpy  and  opaqueness  of  the  gas,  the 
amount  of  material  in  the  layer  of  foam,  the  conductive  heat  flux  to 
the  explosive,  the  heat  evolved  by  the  burning  reaction,  the  enthalpy 
of  the  solid  and  liquid  high-explosive,  and  the  rate  at  which  the  high- 
explosive  reacts  as  a  function  of  temperature.  These  data  have  been 
determined  for  Composition  B  and  are  presented  in  this  report.  The 
problem  of  predicting  the  pressure  rise  during  normal  burning  for 
other  explosive  materials  is  basically  one  of  accumulating  these  data. 
While  some  of  the  data  is  in  the  literature,  most  of  it  must  either  be 
estimated  or  determined  by  experimentation.  Fortunately,  this  is 
not  a  major  endeavor. 

Once  the  high-explosive  cracks,  the  pressure  rises  at  such  a 
rapid  rate  that  much  higher  pressures  are  reached  than  the  system 
can  maintain  under  conditions  of  static  loading.  As  a  result  of  the 
various  modes  by  which  an  explosive  may  crack,  the  rate  of  pressure 
rise  following  the  cracking  (or  the  violence  of  the  resulting  explosion) 
is  not  predictable.  As  a  general  rule  the  order  of  the  explosion  should 
be  related  in  some  inverse  manner  with  the  volume  of  the  void  space 
between  the  explosive  and  the  confining  vessel.  This  is  due  to  the 
fact  that  a  decrease  in  the  void  space  would  accentuate  the  rate  of 
pressure  rise  and  enable  the  confining  vessel  to  remain  intact  until 
higher  pressures  are  reached.  Furthermore,  the  order  of  the  explosion 
should  increase  with  the  depth  of  the  explosive.  The  increased  depth 
would  allow  the  cracking  to  expose  greater  surface  areas  to  the  burning 
and  thereby  result  in  more  pronounced  pressure  rises.  This  would 
partially  explain  the  fact  that  thin  layers  of  explosive  have  never 
initiated  an  explosion  in  any  full-scale  test. 

From  these  studies,  it  is  possible  to  catagorize  the  behavior 
of  any  system  by  one  of  the  following  descriptions: 

a.  Complete  burning  of  the  high-explosive  without  any  hazardous 
pressure  building  up, 

b.  Burning  of  the  high-explosive  that  evolves  a  sufficient  buildup  of 
pressure  to  cause  a  pressure  burst  or  plastic  failure  of  the  confine¬ 
ment  vessel  without  initiating  cracking  of  the  high-explosive.  While 
this  type  of  reaction  would  be  a  hazard  to  personnel  in  the  immediate 
vicinity  of  such  a  reaction,  no  extensive  damage  would  be  anticipated. 

c.  Burning  of  the  high-explosive  that  evolves  sufficient  buildup  of 
pressure  to  crack  the  explosive  material  and  thereby  produce  a  violent 
explosion.  In  this  situation,  the  very  rapid  pressure  rise  may  cause 
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a  brittle  failure  of  the  confinement  vessel  with  the  evolution  of  con¬ 
siderable  shrapnel.  The  explosion  may  or  may  not  be  a  detonation. 

The  criterion  which  determines  whether  (b)  or  (c)  occurs  is  the 
size  of  any  holes  in  the  vessel  at  the  time  the  high- explosive  com¬ 
mences  to  burn.  As  noted  previously  from  the  Composition  B  studies, 
small  holes  of  about  1  percent  of  the  burning  area  will  maintain  the 
pressure  at  near  atmospheric  levels.  With  minor  adjustments,  these 
data  are  also  applicable  to  actual  systems  containing  Cyclotol-type 
explosives.  The  critical  hole  size  will  increase  with  the  proximity 
of  the  hole  to  the  surface  of  the  burning  explosive.  This  is  due  to  the 
lower  density  of  the  gas  caused  by  the  higher  gas  temperatures.  Vent 
holes  should  be  sufficiently  well  dispersed  about  the  high- explosive 
system  to  prevent  any  appreciable  pressure  buildup  as  a  result  of  the 
internal  resistance  presented  to  gas  flow.  Determination  of  the  hole 
sizes  requires  a  knowledge  of  the  parameters  noted  previously  and 
is  a  relatively  simple  procedure. 
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APPENDIX 

Table  A-l  is  a  compilation  of  data  pertaining  to  high-explosives. 
These  data  were  obtained  from  a  literature  study  during  the  course  of 
this  work. 


Explosive 

Material 


Density 
(lb/ ft^) 


'g 


(°F) 

(B/lb) 

Baratol 

(67/33) 

159(a) 

159(b) 

177(j) 

5.  04(J) 

Composition 

B 

101(c) 

103(b) 

1 72 -176(b) 
169-176(c) 
172-176(j) 

14.  4(b) 
14.4(J) 

Cyclotol 

(75/25) 

107(b) 

(cast) 

172-t76(j) 

9.  0(b) 

9.  0(j) 

Cyclotol 

(70/30) 

n 

■ 

Cyclonite 

(RDX) 

1 

114(b) 

399(b) 

HMX 

H 

529(J) 

PETN 

111(b) 

(crystal) 

285(b) 

TNT 

103(b) 

(crystal) 

178(b) 

177.  6(j) 

40.  2(b) 
37.  8(j) 
(calc) 

(a)  ORD  BB.  THI-09D  Pic; 
S^trui?!  Feldman,  Ait.-j: 

(b)  ORDP  20-  177  "Propc  rt 
Military  Interest”  Mav 


TABL 


TABLE  At  PHYSICAL  AND  THERMAL  PROPERTIES  OF  EXPLOSIVES 


.  c,  (B/lb  °F'  and  Thermal  Conductivity,  k,  (B/hr  ft  °F) 
at  Specified  Temperatures  (°F) 

Heat  of 
Explosion 
(B/lb) 

Heat  of 
Formation 
(B/lb) 

■■MM 

> 

32  to  79 

79  to  140 

167 

185 

194 

212 

Cal/Mole 

Temperature 
Ranee  (°F). 

7(j) 

0. 162(a) 

0, 180(j) 

0.280(j) 

0.  213(j) 

0.  201(j) 

0. 171(j) 

0.  193(a) 

0.  248(a) 

0(j) 

0.  254(j) 

0.  305(j) 

0.  376(j) 

0.  354(j) 

0.  341<j) 

0.  312 (J) 

2232(b) 

j 

5(j) 

0.  254(j) 

0.296(j) 

0.  352(J) 

0.  325(j) 

0.  332(j) 

0.  351(J) 

2205(b) 

(calc) 

47,  500(h) 

43,  600(d) 

311.956(d)  1 

8(a) 

5(a) 

0.259(d) 

0.  230(a) 

n.  283(a) 

2304(b) 

0.129(a) 

0.139(a) 

0.  298(b) 

0.  331 

0.  331(b) 
(104°F) 

0.  292(a) 

0.  384(b) 
(176°F) 

0.427(a) 

0.406(b) 

-173(b) 

57,  300(d) 

43.400(d) 

34.400(d) 

47,  500(b) 

338.392(d) 
460-530(d) 
428.500(d) 
415-570  (b) 

0.  288(a) 

0.  288(a) 

1 

3(a) 

0.  262(a) 

0.  262(a) 

0.  282(a) 

2441(b) 

-109(b) 

52.700(b) 

520-5O7(b) 

- 

0.  202(a) 

0.  209(a) 

0.  260(b) 

2493(b) 

689(b) 

50.900(b) 

52,  300(b) 

47,  000(b) 

226-248(b) 

279-314(b) 

322.451(b) 

:<j> 

1(b) 

0.  254(j) 
0.259(a) 

C.  264(d) 

0.  315(j) 

0.  294(a) 

0.  354(j) 

0.  513(j) 

_ 1 

0.478(j) 

0.  275(j) 

1944(b) 

141(b) 

43.400(b) 

34.400(b) 

460-531(b) 

527-590(b) 

0.117(a) 

0,  121(a) 

^  I 

(c)  APG«  OKD  School  APG,  Md,  "Fundamentals  of 
Explosives"  Jan,  1953 

iiio  tviiiianis  and  Wilkes  Co.,  Baltimore  2,  Md. 
"Eighth  Symposium  (International)  on  Combustion"  1962 


(h)  Robertson  A.  J.  B. ,  Trans. 

(j)  R.  Vellcky,  C.  Lenchitz  ant 
"Enthalpy  Change,  Heat  of 
Specific  Heat  of  Basic  Expl 
PATR  No.  2504.  Jan.  1959 


ERMAL  PROPERTIES  OF  EXPLOSIVES 


Vity,  k,  (B/hr  ft  °F) 

Heat  of 
Explosion 
(B/lb) 

Heat  of 
Fc  rmation 
(B/lb) 

ACTIVATION  ENERGY 
<E> 

FREQUENCY  FACTOR  1 
(Z) 

Cal /Mole 

Temperature 
Ranee  <°F) 

1/Sec. 

Temperature 

Rannn 

■1 

■ff mm 

194 

212 

c;280(J> 

0.213(j) 

0.201(j) 

0. 171(J) 

376(j) 

0.  354(j) 

0.  341(J) 

0.  3120) 

2232(b) 

352(j) 

0.  325(j) 

0.  332(j) 

0.  351(j) 

2205(b) 

(calc) 

47,  500(h) 

43,  600(d) 

3U-356(d) 

10)8-  5(h) 

292(a) 

U.  384ib) 
(176°F) 

-  - - 

2304(b) 

2441(b) 

-173(b) 

0.427(a) 

0.  406(b) 

57, 200(d) 

43.400(d) 

34, 400(d) 

47,  500(b) 

3  38  -  392(d) 

460-530(d) 

428-500(d) 
4.15-570  (b) 

1012*19<d) 

1011,4(dl 

1018,5(b) 

460-560(d) 

428.500(d) 

415.570(b) 

288(a) 

282(a) 

-109(b) 

52.700(b) 

520-597(b) 

1019,  7(b) 

520.597(b) 

--  ■ 

— 

2493(b) 

689(b) 

50.900(b) 

52,  300(b) 

47.000(b) 

226-248(b) 

2?9-314(b) 

322.451(b) 

10^0,  b(b) 
io23’  1(b, 

1019*  8(b) 

226-248(b) 

279-314(b) 

322-451(b) 

I54(j ) 

0.  513(j) 

0.478(j) 

0.275(j) 

1944(b) 

141(b) 

43.400(b) 

34,  400(b) 

1 

460-531(b) 

527-590(b) 

1012'2(b) 

1011'  4(b) 

460-531(b) 

527-590(b) 

1 

U>  School  AI’G,  Mu,  "Fundamental!.  of  (h)  Robertson  A.  J.  B. ,  Trans.  Far.  Soc.  45,  85  (1949) 

a"  Jon,  190  3 

(j)  R.  Velicky,  C.  Lenchitz  and  W,  Beach 
ams  ond  Wilkes  Co.,  Baltimore  2,  Md.  "Enthalpy  Change,  Heat  of  Fusion,  and 

ympos ium  (International)  on  Combustion"  1962  Specific  Heat  of  Basic  Explosives" 
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